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ABSTRACT 

Using observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate 
the physical conditions of star-forming regions in z ^ 2.3 galaxies, specifically the electron density 
and ionization state. From measurements of the [O il]AA3726,3729 and [S ii]AA6716,6731 doublets, 
we find a median electron density of ~ 250 cm“^ at z ~ 2.3, an increase of an order of magnitude 
compared to measurements of galaxies at z ^ 0. While z ~ 2.3 galaxies are offset towards significantly 
higher O32 values relative to local galaxies at fixed stellar mass, we find that the high-redshift sample 
follows a similar distribution to the low-metallicity tail of the local distribution in the O32 vs. R23 and 
03N2 diagrams. Based on these results, we propose that z ~ 2.3 star-forming galaxies have the same 
ionization parameter as local galaxies at fixed metallicity. In combination with simple photoionization 
models, the position of local and z ^ 2.3 galaxies in excitation diagrams suggests that there is no 
significant change in the hardness of the ionizing spectrum at fixed metallicity from z ~ 0 to z ^ 2.3. 

We find that z ~ 2.3 galaxies show no offset compared to low-metallicity local galaxies in emission 
line ratio diagrams involving only lines of hydrogen, oxygen, and sulfur, but show a systematic offset 
in diagrams involving [N il]A6584. We conclude that the offset of z ~ 2.3 galaxies from the local 
star-forming sequence in the [N il] BPT diagram is primarily driven by elevated N/0 at fixed 0/H 
compared to local galaxies. These results suggest that the local gas-phase and stellar metallicity sets 
the ionization state of star-forming regions at z ~ 0 and z ~ 2. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: high-redshift 


1. INTRODUCTION 


The star- formation rate density of the universe peaked 


star-forming regions during this epoch is essential for a 
complete description of the evolution of the stellar and 
gaseous content of galaxies. One method of probing these 
conditions is observing recently formed stars embedded 
in ionized gas clouds. Rest-frame optical emission lines 
from HII regions encode a significant amount of informa¬ 
tion about the ionizing source and the physical conditions 
of the ionized gas. A small set of physical properties ap¬ 
pears to determine the strong optical emission line fluxes 
produced in HII regions, including the chemical abun¬ 
dance, the shape and normalization of the ionizing spec- 
trum, the ionization state of the gas, and the gas density 
dPopita fc Ev ans 1986 Kewley & Dopita 2002 iDopita 
et ai.| |2006a|b ). 
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Over the past decade, a body of work has been pro¬ 
duced showing that galaxies at z ~ 1 — 2 display certain 


at z ~ 2 

Hopkins & Beacom 

2006 

Madau & Dick- 


inson 

2014 

. Understanding the physical conditions in 


galaxies (e.g., Shapley et al. 2005 

|Erb et al.| 20061 Liu 

et al.|20U81|Hamline et al.|2U091|Ho 

iden et ai.|2[jl41|iNaka- 

iima & Ouchi||2014| |Steidel et al.||2014| |Shapley et al. 


2015|). These observations suggest that at least some 


of the physical conditions that influence emission line 
production in star-forming regions must be different in 
high-redshift galaxies. It is well-established that galax- 


at fixed stellar mass ( 

Erb et al. 2006 Maiolino et al. 

2008 

Belli et al.| 

2013 

IHenry et al.||2013t IStott et al. 

2013 

Cullen et a 


141 jSteidel et al. 2014 Troncoso 

et al. 

2014 Sanders et ai.| 20151. The electron density 


production, but are less well-constrained and have fewer 
measurements at z > 1. 

Estimates of the electron density at z ^ 2 have sug¬ 
gested that densities in high-redshift star-forming regions 
are sign ificantly higher than what is typically observed 


locally (Hainline et al. 

200^ 

Lehnert et al. 

2009 

iBian 

et al.|20l0 bhirazi et aJ 

.112014 

I. However, the small and 

heterogeneous galaxy samples 

used in these studies leave 


the electron density of typical star-forming regions at 
z ^ 2 poorly constrained. Electron density estimates 
in a large, representative sample of z ~ 2 galaxies are 
needed to obtain robust constraints on the typical elec¬ 
tron densities at high redshift. 

The typical ionization state of gas in z ~ 2 
star-forming regions appears to differ from that 
of local star-forming regions. Galaxies at z > 

I display higher values of [O iii]A5007/H/3 and 
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[O iii]AA4959,5007/[O ii]AA3726,3729 than are typical 
of local galaxies, leading to the conclusion that high- 
redshift galaxies have hig her ionization parameters than 
local galaxies, on average dBrinchmann et al.|2008 |Hain- 
line et aL]|2009[ IMolden et al.||2U14[ [INakajima fc Uuchi, 
2014 Shirazi et al.||2014 ). A systematic investigation of 
the variation ot the ionization state with other galaxy 
properties at z ~ 2 is necessary to uncover the cause of 
these elevated ionization parameters. 

In this paper, we investigate the physical proper¬ 
ties of star-forming regions at z ~ 2.3, specifically 
the electron density and ionization state, using a large, 
systematically-selected sample from the MOSFIRE Deep 
Evolution Field (MOSDEF) survey. Until recently, sam¬ 
ples of rest-frame optical spectra of z ~ 2 galaxies were 
small, heterogeneous, and usually included only a subset 
of the strongest nebular emission lines. These samples 
reflected the difficulty of obtaining spectra of faint galax¬ 
ies in the near-infrared, typically using long-slit spectro¬ 
graphs covering one near-infrared band at a time. With 
the development of sensitive near-infrared detectors and 
multi-object near-infrared spectrographs on 8-10 m class 
telescopes, large samples of z ~ 2 galaxies with rest- 
frame optical emission line measurements across multiple 
near-infrared bands are being assemb led for the first time 
(Steidel et al. 2014 Kriek et al. 20151. These spectra con- 
tain a wealth of diagnostic information that probes the 
star-formation rate (SFR), dust attenuation, gas density, 
ionization state, chemical enrichm ent, and more fo r z ~ 2 


galaxies. Building on the work of Shapley et al. (2015), 
we utilize the full range of strong optical emission lines 
to investigate the physical properties of HII regions using 
various diagnostic line ratios. In combination with rich 
datasets at lower redshifts, such measurements make it 
possible to understand how conditions in star-forming 
regions have evolved over the past ~ 10 Gyr of cosmic 
history. 

Changes in the physical conditions of star-forming gas 
are thought to be the cause of the offset of z > 1 galax¬ 
ies from the local sequence of star-forming galaxies in 
the [O iii]A5007/H/3 vs. [N ii]A6584/Ha excitation dia¬ 


gram (iSJrapley et al.|2005l|Liu et al.|2008 
2009| ISteidei et ai.|[2U14| Shapley et al. 
been proposed th at tiiis offset is cause d' 


Hainline et al. 


^). It has 

__ by higher gas 

density/pressure dKewl ey et al.||2013a|), sy stematically 
harder ionizing spec tra (fsteidel et al.||2014 b. higher ion¬ 
ization parameters (Brinchm ann et al. 2008), or an en¬ 
hancement in t he N/U ratio (|iVlasters et al.||2014| |Shap- 
ley et al. 2015) in high-redshitt galaxies m comparison 
to what is typically observed in the local universe. The 
offset could be caused by a combination of some or all 
of these parameters. By characterizing the differences 
in the density and ionization state of z ^ 2 and z ~ 0 
star-forming galaxies, as we do here, we can gain a bet¬ 
ter understanding of which parameters drive the offset 
in the [O iii]A5007/H^ vs. [N ii]A6584/Ha excitation 
diagram, and the relative importance of each. Determin¬ 
ing the evolution of these properties with redshift also 
has implications for the applicability of local metallicity 
calibrations at z ~ 1 — 2. 

This paper is structured as follows. In Section we 
briefly describe the MOSDEF survey, along with the ob¬ 
servations, reduction, and measurements. We estimate 
the typical electron density in z ~ 2.3 star-forming re¬ 


gions and characterize the evolution of density with red- 
shift in Section In Section we investigate the ion¬ 
ization state of z ~ 2.3 galaxies and its dependence on 
global galaxy properties and metallicity indicators. We 
propose a scenario in which galaxies at z ~ 2.3 have 
the same ionization parameter as galaxies at z ^ 0 
at fixed metallicity. In Section we provide evidence 
supporting our proposed scenario and discuss the impli¬ 
cations for the interpretation of observed emission-line 
ratios of z ^ 2.3 galaxies, including the offset in the 
[O iii]A5007/H/3 vs. [N ii]A6584/Ha diagram. We con¬ 
clude by summarizing our main results in Section]^ 

We adopt the following shorthand abbreviations to re¬ 
fer to commonly used emission line ratios: 

O32 = [O iii]AA4959,5007/[0 ii]AA3726,3729 (1) 

R23 = ([O iii]AA4959,5007-f [O ii]AA3726,3729)/H/3 

( 2 ) 

03N2 = ([O iii]A5007/H/?)/([N ii]A6584/Ha) (3) 

N2 = [N ii]A6584/Ha (4) 

Throughout this paper, the term “metallicity” is used 
synonymously with gas-phase oxygen abundance (0/H) 
unless otherwise mentioned. We adopt a A-CDM cos¬ 
mology with Hq = 70 km s“^ Mpc“^, Vim = 0.3, and 
Da = 0.7. 


2. OBSERVATIONS 

We use data taken during the first two years (2012B- 
2014A) of the MOSFIRE Deep Evolution Field survey. 
We briefly describe the MOSDEF survey, observations, 
reduction, and derived quantities here. Full technical 
details of the survey strategy, observations, reduction 
pipeline, and sam ple characteristics can be found in 
Kriek et al. (2015). We additio nally use data fro m the 
Sloan Digital Sky ^Survey (SPSS; [York et al.||2000|) Data 
Release 7 (DR7; |Abazajian et al.||2009l) catalog to se¬ 
lect local comparison samples for studying evolution with 
redshift. Emission-line measurements and galaxy prop¬ 
erties are taken from the MPA-JHU catalog of measure¬ 
ments for SDSS DR70 


2.1. The MOSDEF Survey 

The MOSDEF survey is an ongoing multi-year project 
in which we are obtaining rest-frame optical spectra of 
galaxies at z~ 1.4 — 3.8 with the goal of transforming 
the understanding of the gaseous, stellar, dust, and black 
hole content of galaxies at that epoch in cosmic history. 
This project utilizes the Multi-Obj ect Spectrometer For 
Infra-Red Exploration (MOSFIRE; McLean et al.||2012 1 
on the 10 m Keck I telescope. Potential for scientific gam 
from the MOSDEF dataset is maximized by targeting 
objects in the AEGIS, GOSMOS, and GOODS-N extra- 
galactic fields with extensive multi-wavelength ancillary 
data. These data include Hubble Space Telescope (HST) 
imaging from the Gosmic Assembly Near- infrared Deep 
Extragalactic Legacy Survey ( CANDELS; |Grogin et al.| 
2011 Koekemoer et al. 20111) and grism spec troscopy 
from the 3D-11S4’ survey ^Brammer et al.|2012a ), as well 
as observations from Chandra, Svitzer, Herschel, VLA, 


Available at http://www.inpa-garching.inpg.de/SDSS/DR7/ 
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and ground-based observatories in the optical and near- 
infrared. 

In the MOSDEF survey, we target galaxies in the three 
redshift windows 1.37 < z < 1.70, 2.09 < z < 2.61, and 
2.95 < z < 3.80, where the redshift ranges are selected 
such that strong optical emission-line features fall within 
windows of atmospheric transmission in the Y, J, H, or K 
near-infrared bands. This targeting strategy leads to cov¬ 
erage of [Oii]AA3726,3729, H/3, and [O iii]AA4959,5007 
for all three redshift bins, as well as [N ii]AA6548,6584, 
Ha, and [S ii]AA6716,6731 for the z ~ 1.5 and z ~ 2.3 
bins. These strong optical emission lines are powerful 
diagnostics for determining physical properties of galax¬ 


ies including dust cont ent (Reddy et al. 
formation r ate (SFR) jShivaei etjF ' 


abundance (Sanders et aT 2015), gas density, ionization 


in 


2015), star- 


0 ), cnemical 


state, and black hole activity (|Coil et al.||2015|). Targets 
are selected fro m the 3D-HST phot ometric and spectro¬ 
scopic catalogs ( Skelton et al.||2014 ) based on their rest- 
frame optical (observed H-bandj magnitudes and red- 
shifts (grism or spectroscopic redshifts, if available, and 
photometric redshifts otherwise). Galaxies are targeted 
down to HST /WFC3 F160W AB magnitudes of 24.0 at 
z ^ 1.5, 24.5 at z ^ 2.3, and 25.0 at z ^ 3.4. Tar¬ 
gets with pre-existing spectroscopic redshifts, including 
grism redshifts or ground-based spectroscopic redshifts, 
are given higher priority, as are gal axies with brighte r 
F160W magnitudes. As described in Kriek et al. (2015), 
the MOSDEF survey will obtain spectra for ~ 1500 
galaxies when complete, with ~ 750 at z ^ 2.3, and 
^ 400 each at z ~ 1.5 and z ~ 3.4. 

2.2. Observations and Reduction 

We use data from the first two observing seasons of 
the MOSDEF survey. Observations were taken on ten 
observing runs from 2012 December to 2014 May, dur¬ 
ing which 21 MOSFIRE masks were observed. The first 
observing run, taking place in 2012 December, was a pi¬ 
lot program during which we observed one mask each 
in the GOODS-S and UDS GANDELS fields due to the 
limited visibility of the primary target fields. This work 
focuses on the z ^ 2.3 redshift bin, and, accordingly, we 
only describe the observations for this redshift interval. 
MOSDEF targets at z ^ 2.3 are observed in the J, H, 
and K near-infrared bands, with [Oil]AA3726,3729 in J, 
H/3 and [O iii]AA4959,5007 in H, and [N ii]AA6548,6584, 
Ha, and [S ii]AA6716,6731 in K. Observed MOSFIRE 
masks each contain ~ 30 slits with widths of 0'.'7, yield- 


servations of B8-A1 V standard stars matched to the 
typical air mass of science observations. Flux calibra¬ 
tion was performed by requiring the flux density of the 
reference star on a mask to match its cataloged broad¬ 
band photometry. For each slit, a two-dimensional error 
spectrum was produced accounting for Poisson count¬ 
ing uncertainty for the observed intensity per pixel and 
read noise. From the two-dimensional science and er¬ 
ror spectra, one-dimensional spectra were produced us¬ 
ing the optimal extraction technique. Spectra for any 
detected objects serendipitously falling on the slit were 
also extracted (Freeman et al., in prep.). The final flux 
calibration was achieved by applying a slit-loss correction 
term to the extracted science spectra on an individual ba¬ 
sis. The fraction of light from an object falling outside of 
the slit was estimated using a two-dimensional Gaussian 
fit to the F160W image of a galaxy convolved with the 
seeing estimate for each mask and filter. The flux calibra¬ 
tion was checked by verifying that objects with detected 
continuum had flux densities consistent with broadband 
photometry. 

2.3. Measurements and Derived Quantities 

Emission line fluxes were measured by fitting Gaussian 
profiles to emission lines in the one-dimensional spec¬ 
tra, while the uncertainty on the emission line flux was 
based on the 68th percentile width of the distribution 
of measured fluxes obtained by perturbing the spectrum 
according to the error spectrum and refi tting the emis¬ 


sion line 1,000 times (Reddy et al. 2015). All emission 
lines were fit with a single Gaussian profile except for 
the cases of [N ii]AA6548,6584 and Ha, fit simultane¬ 
ously with a triple Gaussian, and [O ii]AA3726,3729, fit 
simultaneously w ith a double Gaussian and described in 
detail in Section |3.2| The highest signal-to-noise emis¬ 
sion line of each object, typically Ha or [O iii]A5007, was 
used to constrain the centroid and width of the other 
emission lines and measure the redshift. Stellar masses 
were estimate d by utilizing the stel lar population synthe¬ 


sis models of Gonroy et al. (2009) with the SED-fitting 
code FAST ( Kriek et af .|2009 ) using the measured MOS¬ 


DEF spectroscopic redshifts and broadban d photometric 


catalo gs assembled by the 3D-HST team (Skelton et ah, 
[20] ^ spannin g observed optical to mid-infrared. I'he 


Galzet ti et al. (2000) attenuation curve and a Chabrieii 
(20U3| IMh' were assumed. Uncertainties on the stel- 


lar masses were estimated using a Monte Garlo method 
where the input photometry was perturbed according to 


and K bands, respectively. One slit on each mask was 

2015). SFRs were estimated from dust-corrected Ha fu- 

placed on a reference star used in the reduction. Masks 

minosities using the 

Kennicutt 

(1998) calibration con- 

were typically observed using an ABA’B’ dither pattern, 

verted to a 

Chabrier 

120931 IM 

F (Shivaei et al. 

2015 

). 


Balmer line fluxes were first corrected tor steffar Balmer 


120 seconds in J and H, and 180 seconds in K. The total 
exposure time per filter per mask was typically 2 hours. 

The raw data were reduced using a custom IDL 
pipeline de veloped by the MOS DEF team and described 
in detail in 
were flatfie' 


absorption using the best-fit SEDs (Reddy et al. 


2015) 


Kriek et al. (2015). The raw science frames 
ded and sky subtracted, cosmic rays were 
identified and masked, and the two-dimensional spectra 
were rectified. Individual exposures were combined and 
the resulting spectrum was flux calibrated. Shape cor¬ 
rection due to varying spectral response with wavelength 
and telluric absorption features was achieved using ob- 


and the dust correction was estimated usi ng Ha/H/3 as¬ 
sumin g an intrinsic ratio of 2.86 and the [Cardelli et ah] 
(19891 extinction curve. 

Emission line fluxes are dust-corrected in line ratios in¬ 
volving lines significantly separated in wavelength (O32 
and R23), but are not dust-corrected in line ratios fea¬ 
turing lines with a close wavelength spacing (03N2, N2, 
[O iil]/H/3, and [S ii]/Ha). Uncertainties on emission 
line ratios are estimated using a monte carlo method by 
perturbing each individual line flux according to its un- 
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certainty, recalculating the emission line ratio, and re¬ 
peating the preceding steps 10,000 times to build up a 
distribution of perturbed line ratios. The uncertainty 
on the line ratio is determined from the 68 th percentile 
width of this distribution. Error estimates on O 32 and 
R23 values include uncertainty in the Balmer decrement 
by recalculating the dust correction for each realization. 
After removing those objects i dentified as AGN based 


on their X-ray or IR properties (Coil et al. 2015), there 
are 225 star-forming galaxies conhrmed to be in the 
z ~ 2.3 redshift interval spanning a stellar mass range 
of 10 ®'®^ — 10 ^^ ®"^ Mq with a median stellar mass of 
109.99 subset of these galaxies with Ha and 

H/? detected (67%) has SFRs spanning 1.61 — 323 MQ/yr 
with a median SFR of 21.6 MQ/yr, and spans a range in 
stellar mass of 10 ®- ^^ — 10 ^^-^^ Mq with a median stellar 


mass of 10^®'® Mq. Shivaei et al. (2015) have shown that 
the SFR-M* relation of MUSUFF z ~ 2.3 star-forming 
galaxies is similar to what is observed in other studies 
that employ different SFR indicators (see their Figure 8 ). 
Accordingly, the z ~ 2.3 star-forming galaxy population 
recovered by the MOSDEF survey appears to be repre- 


galaxies at these stellar masses and this redshift (Reddy 

et al .||2012 

Whitaker et al.||2014 

Shivaei et al.||2015p. 


3. ELECTRON DENSITIES 

The electron density of star-forming regions affects the 
fluxes of collisionally excited lines and is thus an impor¬ 
tant quantity to measure as an input parameter to pho¬ 
toionization models^ The electron density can be esti¬ 
mated using the line fluxes of components of a doublet 
of a single species in which the two levels of the doublet 
have different collision strengths and radia tive transition 
probabilities ( Osterbrock fc Ferland|[2006 ). The flux ob¬ 
served in eacn component of the doublet is dependent 
on the relative population in each energy level, which is 
sensitive to the density of electrons available for colli- 
sional excitation and de-excitation. Rest-frame optical 
spectra provide access to two strong emission line dou¬ 
blets useful for estimating the electron density, namely 
[O ii]AA3726,3729 and [S ii]AA6716,6731. We use the ra¬ 
tios of the components of these doublets to estimate den¬ 
sities for local SDSS galaxies and high-redshift galaxies 
from the MOSDEF sample. 


3.1. Methods 

We have written a python script that solves a 5- 
level atom approximation of the O ii and S ii ions for 
the relative populations in the second and third en¬ 
ergy levels. Decays from these two energy levels pro¬ 
duce the two components of the [O ii]AA3726,3729 and 
[S ii]AA6716,6731 doublets. We set up a detailed bal¬ 
ance of transitions into and out of each of the five energy 
levels via radiative decay and collisional excitation and 
de-excitation, assuming the system is in thermal equilib¬ 
rium. The detailed balance provides a system of equa¬ 
tions that can be solved for the relative populations at a 
given density. Given 712 and 773 , the relative populations 

^ The models described in Section!^ actually use the hydrogen 
gas density as an input parameter. In HII regions, the hydrogen 
gas is fully ionized such that the electron density is a good proxy 
for the hydrogen gas density. 


in the second and third energy levels, the ratio of the line 
fluxes is given by 


P3 _ ^ 773^31 

F2 F21772A21 n 2 A 2 i 


( 5 ) 


where Fi is the emission line intensity of decay from the 
7 th level to the ground state, Fn is the energy difference 
between the fth level and the ground state, and An is 
the transition probability of the fth level to the ground 
state. 

Calculating the proper emission line ratio correspond¬ 
ing to a given density requires accurate knowledge of the 
transition probabilities and collision strengths of each 
transition between the five energy levels. Recent inves¬ 
tigations (Nicholls et al. 2013 Dopita et al. 2013) have 
suggested advantages in using the most up-to-date col¬ 
lision strength and transition probability atomic data 
instead of o utdated values include d in the IRAF rou¬ 
tine temden ( Shaw fc Dufour|[i994 1. Motivated by these 
studies, we adopt the eff ective collision strengths f rom 
Tayal (2007) for O ii andjTayal & Zatsarinny (2010) for 
S II, while the transition probabilities tor b oth species 
are taken from the NIST MCHF database (Fischer & 
Tachiev 2014). Using other estimates of the collision 
strengths and transition probabilities can change the cal¬ 
culated electron densities by up to ^ 30%. We verified 
that our script can exactly match IRAF temden when 
using the same atomic data. 

The effective collision strengths have some tempera¬ 
ture dependence and have been calculated over a range 
of electron temperatures from 2000 K to 100,000 K for 
O II and 5000 K to 100,000 K for S ii. We adopt the effec¬ 
tive collision strengths calculated with an electron tem¬ 
perature of 10,000 K, a representative equilibrium tem¬ 
perature of HII regions that are neither metal-rich nor 
metal-poor. We note that the temperature dependence 
is not negligible. Assuming an electron temperature of 
7,000 K yields electron densities that are ~ 15 — 20% 
lower at a fixed line ratio, while an electron temperature 
of 15,000 K gives densities that are higher by the same 
amount. Because we assume a fixed electron tempera¬ 
ture of 10,000 K, we are likely overestimating the elec¬ 
tron density in metal-rich galaxies, and underestimating 
the density in galaxies that are metal-poor. However, 
the uncertainty introduced by this assumption is smaller 
than the typical measurement uncertainty for individual 
z ^ 2.3 galaxies. 

We calculate the line ratios [O ii]A3729/A3726 and 
[S ii]A6716/A6731 over a range of electron densities of 
log = 0 to 5 in 0.01 dex steps. The result is well 

fit by a function of the form 


Ririe) = 


( 6 ) 


where ii]A3729/A3726 or [S ii]A6716/A6731 is the 
line flux ratio. The best-fit parameters using up-to-date 
atomic data are shown in Tablej^for [O il] and [S il]. Fig¬ 
ure [ 1 ] shows this diagnostic relation for the new atomic 
data (black) and the relation from IRAF temden (red) 
for both [O ii] (solid) and [S ii] (dashed). It can be 
seen in Figure that the line ratio asymptotically ap¬ 
proaches a theoretical maximum value in the low-density 
limit and a theoretical minimum value in the high-density 
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Table 1 

Coefficients and limiting line ratios for [O ii] and [S ii] in 
equations [^and[^ 


R 

a 

b 

c 

p a 

-‘Tmin 

itmax® 

[O ii]A3729/A3726 
[S ii]A6716A6731 

0.3771 

0.4315 

2,468 

2,107 

638.4 

627.1 

0.3839 

0.4375 

1.4558 

1.4484 


^Theoretical minimum line ratio calculated in the high-density 
limit of 100,000 cm“^ 

'^Theoretical maximum line ratio calculated in the low-density 
limit of 1 cm“^ 



Figure 1. R vs. rig curves (equation from IRAF temden 
(red) and our five-level atom python script using new atomic data 
(black), where R =[0 ii]A3729/A3726 (solid) or [S ii]A6716/A6731 
(dashed). 


limit. We calculate the theoretical maximum and mini¬ 
mum line ratios at densities of 1 cm“^ and 100,000 cm“^, 
respectively, and show these values in Table Inverting 
equation yields the density as a function of the line 
ratio 


ne{R) 


cR — ab 
a — R 


(7) 


which we use to calculate electron densities. Uncertain¬ 
ties on individual density measurements are estimated 
by converting the upper and lower 68 th percentile un¬ 
certainties on the line ratio into electron densities, where 
the upper (lower) uncertainty in line ratio corresponds 
to the lower (upper) uncertainty in density. 


3.2. Sample 

Eor objects with electron densities in the range ^ 
100 — 1,000 cm“^, the line ratio will be fairly close to 
unity, as can be seen in Figure In this regime, rela¬ 
tively small changes in the line nux can result in large 
changes to the inferred electron density. Therefore, it is 
imperative that the doublets used to infer electron den¬ 
sities be free of any contamination from skylines or poor 
line profile fitting. The [S ii] doublet is well separated 
and the two components are fit separately with single 


Gaussian prohles. The [O il] doublet is well-resolved but 
the two components are blended and must be fit with two 
Gaussian prohles simultaneously. When htting the [O il] 
doublet, we constrain the separation of the centroids of 
the two components to be within 0.5 A of the nominal 
separation of 2.78 A in the rest frame. We also require 
the widths of the two components to match each other 
exactly and to be no more than 10 % larger than the ve¬ 
locity width inferred from the highest signal-to-noise line 
in the object’s spectrum, typically Ha or [O iii]A5007. 
This method yields robust hts to [O ii] doublets. 

We select MOSDEF galaxies in the redshift range 
2.0 < z < 2.6 with S/N>3 in [O ii]AA3726,3729 or 
[S ii]AA6716,6731 which have not been flag ged as AGN 


based on their IR and X-ray properties (Coil et al. 
This gives a sample of 97 [O ii] doublets and 


2015 

36 [S ii] doublets at z ~ 2.3. We visually inspected 
each of these doublets and removed those with sig¬ 
nificant skyline contamination or spurious detections. 
One additional [S ii] object with a very high value of 
log ([NiiA6584/Ha) = —0.11 indicating a probable AGN 
was also removed, giving a final density sample of 43 
[O ii] doublets and 26 [S ii] doublets from 61 different 
targets at (z) = 2.24 ± 0.12. This sample spans a range 
of stellar mass from 10 ®'®^ — 10 ^^’^^ Mq with a median 
stellar mass of 10 ^°-^° M©. Of the 61 galaxies, 8 (13%) do 
not have measured SFRs, of which 3 galaxies do not have 
wavelength coverage of H/3, 2 galaxies do not have wave¬ 
length coverage of Ha, and the remaining 3 galaxies show 
significant skyline contamination of the H/3 line. The 
other 53 objects have SFRs spanning 4.65 — 228 M 0 /yr 
with a median SFR of 29.7 M 0 /yr. The density sample 
has only slightly higher median M* and SFR than the 
parent sam ple o f MOSDEF z ^ 2.3 star-forming galax¬ 
ies (Section 2.3) and is still representative of the fairly 
massive z ^ 2.3 star-forming galaxy population. The 
8 galaxies without SFR measurements are included in 
all parts of the analysis for which SFR is not required. 
Examples of [O ii] and [S ii] doublets from six different 
objects in the MOSDEF z ^ 2.3 sample are shown in 
Figure]^ along with fits and inferred line ratios. 

We use a local comparison sample from SDSS to inves¬ 
tigate evolution in the typical densities of star-forming 
regions. We select galaxies from SDSS in the redshift 
range 0.04 < z < 0.1 to attain a sample that is rel¬ 
atively free of aperture effects and limited to the local 
universe. We require galaxies to have S/N>3 in H/3, 
[O iii]A5007, Ha, [N ii]A6585, [S ii]A6716, and [S ii]A6731. 
Detections in the first four of these lines are requ ired to 
reject AGN using the empirical demarcation of |Kauff-| 


mann et al. (2003). We do not require detection of the 
[U II] doublet, ihe SDSS spectra have a spectral reso¬ 
lution of ^ 2000 , corresponding to a resolution element 
of ^ 1.9 A at 3727 A. This resolution is insufficient to 
properly resolve and sample the [O ii]AA 3726,3729 dou¬ 
blet separated by 2.78 A (see Section 3.3). Therefore, we 
only use the [S ii] doublet to probe the electron density 
in the local comparison sample. The local comparison 
sample contains 99,291 galaxies with (z) = 0.0678. 


3.3. The consistency of [O ii] and [S ii] electron 
densities 
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Figure 2. [O ll]AA3726,3729 doublets (left column) and 
[S ll]AA6716,6731 doublets (right column) from six different ob¬ 
jects over a range of line ratios and densities. The black line shows 
the continuum subtracted spectrum for each object. The light gray 
band indicates the error spectrum for each object, while the blue 
and red lines show the Gaussian profile fits to the blueward and 
redward component of each doublet, respectively. The green line 
shows the total [O ll] profile. In each panel, the spectrum has been 
normalized so that the blue component has a peak height of unity. 
The line ratio R, density (ue) in cm“®, and redshift is displayed 
for each object, with R =[0 ll]A3729/A3726 in the left column and 
R =[S ll]A6716/A6731 in the right column. 


Since our high-redshift electron density sample has 
a mixture of [O il] and [S il] doublet measurements 
while the local comparison sample only has reliable 
[S ii] measurements, it is a useful exercise to evalu¬ 
ate the consistency of densities determined using these 
two ionic species to see if they can be directly com¬ 
pared. To this end, we have assembled a sample of lo¬ 
cal HII regions from the literature with high-resolution, 
high-signal-to-noise spectroscopic observations with suf¬ 
ficient wavelength coverage to span [O ii]AA3726,3729 
and [S ii]AA6716,6731. We performed a literature search 
and identified 32 galactic and extragalactic HII regions 
observed at high spectral resolution (i? ~ 8,000—23, 000) 
with detections of both the [O i i] and [S ii] doublets 


(|Garcia-Ro 

as et al. 2005 2006 12007 

Lopez-Sanchez 

et al.||2007' 

Esteban et al.||2009 

20T3'2D] 

H). 


Electron densities and uncertainties are calculated with 
the same method outlined above using the published line 


fluxes and errors. Densities and uncertainties are pre¬ 
sented in Figure Four of the 32 individual measure¬ 
ments have [S ii]A6716/A6731 ratios that are higher than 
the theoretically allowed maximum in the low-density 
limit, and thus cannot be assigned a density. While these 
four objects have [S ii]AA6716,6731 detected at greater 
than 3 ct indicating they are in the low-density limit, we 
plot them as upper limits (red squares) where the data 
point is plotted at the [S il] density corresponding to the 
lower Icr uncertainty on the [S ii] line flux ratio. Fitting 
a line in logarithmic space yields the relationship 


where the 68th percentile confidence intervals are deter¬ 
mined by perturbing the data points according to their 
uncertainties and refitting. This best-ht line is shown in 
Figure as a solid blue line while the light blue shaded re¬ 
gion shows the 68th percentile confidence region around 
the best-fit line. The relation between electron densities 
determined by [O ii] and [S ii] is completely consistent 
with a one-to-one relation (dashed black line). All four 
objects plotted as upper limits are also consistent with 
a one-to-one relation within the Icr uncertainties. There 
are 8 objects in the z ~ 2.3 density sample that have 
density estimates from both the [O ii] and [S il] doublets, 
but the number of galaxies is too small and measurement 
uncertainties are too large to perform a similar investiga¬ 
tion at z ~ 2.3. We assume that the relationship between 
densities of star-forming regions determined using [S il] 
and [O If] does not change with redshift. Thus, we con¬ 
clude that densities determined from either ionic species 
in the z ^ 2.3 sample can be directly compared with each 
other and with SDSS density measurements from [S il] 
doublets. 

For this test, we specifically selected a sample with 
high S/N and very high spectral resolution so that each 
component of the doublets was well-detected and the 
[O ii] doublet was well resolved. We note that repeat- 


1,000 — 2,000) of local HII regions 1 

Peimbert et al.|2012 

Garcfa-Rojas et al.|2014 

Berg et a 

.||2015|) yields a rela- 


overestimated with respect to [S ii] electron densities. 
This effect is very similar to what is seen in the SDSS 
sample {R ~ 2, 000). In order to have at least two resolu¬ 
tion elements to sample the separation of the [O ii] dou¬ 
blet components, a spectral resolution of AA = 1.39 A 
at A = 3727 A is needed, corresponding to i? ~ 2, 700. 
MOSDEF observations adequately sample the [O ii] dou¬ 
blet at z ^ 2.3 with i? ~ 3, 300 in the J band. 


3.4. Typical electron density at z 2.3 and z ~ 0 

We would like to characterize the typical electron den¬ 
sity in star-forming regions of z ^ 2.3 galaxies and z ~ 0 
galaxies. However, given the shape of the function in 
equation and Figure a fairly symmetric distribu¬ 
tion of line ratios leads to a very asymmetric distribu¬ 
tion in electron densities. Furthermore, the diagnostic 
curve translating from line ratio to electron density is 
insensitive to the electron density at very low and very 
high densities, asymptotically approaching the theoret¬ 
ical maximum and minimum line ratio, respectively, in 
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Figure 3. Comparison of density estimates from the [O ll] and 
[S ll] doublets for a sample of local HII regions with high-S/N, high- 
resolution spectra. Black points denote density measurements for 
individual HII regions. The four red squares show limits plotted 
at the upper Icr uncertainty bound on the [S ll] density for objects 
that have higher [S ll]A6716/A6731 than the maximum theoretically 
allowed value. The black dashed line shows a one-to-one relation¬ 
ship. The blue line and shaded blue region show the best-fit line 
and IfT confidence interval, respectively. Parameters of the best-fit 
line are shown in equation]^ 

those two regimes. Measured line ratios that fall outside 
of the theoretically allowed region due to measurement 
uncertainty can only be assigned limits in the low- or 
high-density extremes. We consider the low-density limit 
to refer to a density below ~ 10 cm“^ and the high- 
density limit to denote density above ^ 10,000 cm“^. 
For these reasons, we perform statistics on the line ra¬ 
tio distributions for each sample rather than the electron 
density distributions and infer typical electron densities 
based on the statistical properties of the line ratio distri¬ 
bution. 

The distributions of line ratios of the z ~ 2.3 star¬ 
forming galaxies are shown in the top and bottom panels 
of Figure for the [O il] and [S il] doublets, respec¬ 
tively. We determine the typical electron density from 
the median line ratio of a given sample. The uncertainty 
on the median is calculated using a bootstrap technique 
in which we randomly resample with replacement, per¬ 
turb the emission line fluxes according to their uncer¬ 
tainties and recalculate the line ratio for each object in 
the new sample, take the median of the new perturbed 
sample, and repeat the preceding steps 1,000 times to 
build up a well-sampled distribution of median values. 
The reported lower and upper uncertainties on the me¬ 
dian are determined to be the 15.8-percentile and 84.2- 
percentile values, respectively, of the cumulative distri¬ 
bution function of the median. At z ~ 2.3, we find a me¬ 
dian [O ii]A3729/A3726 ratio of 1.181 q; 5J corresponding 
to an electron density of 22514^® cm“^. We find a me¬ 
dian [S ii]A6716/A6731 value of 1.13lo;|!)g which gives an 
electron density of 290^169 cm consistent within the 



Figure 4. [O ii]A3729/A3726 (top) and [S ii]A6716/A6731 (bot¬ 
tom) line ratio distributions for 43 and 26 2 : ~ 2.3 star-forming 
galaxies, respectively. In each panel, the dotted black line shows 
the median line ratio (corresponding to an electron density of 
225 cm“^ for [O ll] and 290 cm“^ for [S ll]), while the dashed 
red lines show the minimum and maximum theoretically-allowed 
line ratios from Tabled 



[S II]A6716/A6731 


Figure 5. [S ii]A6716/A6731 line ratio distribution for local star¬ 
forming galaxies from SDSS, with lines as in Figurel^ The median 
line ratio for local galaxies (corresponding to an electron density 
of 26 cm“^) falls near the low-density limit. 
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uncertainties with the density determined using [O ii]. 
We measure a range of individual electron densities from 
the low-density limit to 2,500 cm“^ and hnd that z ~ 2.3 
star-forming regions have a typical electron density of 
^ 250 cm“ . 

The typical density that we infer for high-redshift 
galax ies 250 cm~^) is in excellent agreement with 
what iSteidel et al.| (|2014|) observed by stacking J-band 
observations of 113 galaxies at z ~ 2.3, finding an av¬ 
erage [O il]A3729/A3726 ratio of 1.16, corresponding to 
a density of 2 43 cm~^ using th e atom ic data adopted 


in this paper. Shimakawa et al. (2015) found a median 
electron density of 291 cm“'^ among 14 Ha emitters at 
z ^ 2.5 using the [O ii] doublet. Previous observations 
of individual gravitationally lensed galaxi es at z ^ 2 sug¬ 
gested electron densities of ^ 1000 cm“^ (Hainline et al. 


2009 iBi an et all2010|), somew hat higher than the value 
we inter. [Lehnert et al. (2009) presented electron densi¬ 
ties in the range 400 —1200 cm~^ for 4 galaxies at z ~ 2.3 
in the SINS survey (iForster Schreiber et al.||2009D, also 


higher than our sample median. As noted m Section |3.1 
different choices of atomic data can change the inferrec 
densities by ~ 30% and can lead to differences of this 
magnitude in measured densities reported by different 
authors. In co mparison to these pr evious estimates, ex¬ 
cept for that of|Steidel et al.|(|2014|), our sample is larger 
and selected in a more systematic way, and the galax¬ 
ies in our sample display properties representati ve of the 
SFR-M, relation at z ~ 2, as shown in Section [23] As 
such, the density estimate presented here should hold 
true for a population of typical star-forming galaxies at 
z ~ 2.3 with M* > 10®'5 Mq. 

For the local comparison sample, we find SDSS star¬ 
forming galaxies have a distribution with a median 
[S il]A6716/A6731 ratio of 1.41, shown in Figurej^ corre¬ 
sponding to an electron density of 26 cm“^. The uncer¬ 
tainty on the SDSS median line ratio is less than 0.04% 
due to the large number of galaxies in the sample. We 
find that the typical electron density in star-forming re¬ 
gions increases by a factor of 10 from z^0toz~2.3. 
The local median [S ii]A6716/A6731 ratio is close to the 
theoretical maximum ratio of 1.4484, suggesting that lo¬ 
cal star-forming galaxies typically fall close to the low- 
density limit. From the slope of the function in Figure!^ 
it is apparent that the line ratios of both [O ii] and [S li] 
are almost completely insensitive to the density when 
the electron density is below ~ 10 cm“^, and are only 
mildly sensitive to the density below ~ 50 cm“^. Even 
if we assume a conservative upper limit of 50 cm“^ for 
the typical local density, we still observe a significant in¬ 
crease in electron density from z~0toz^2.3. 

The measurements above suggest that z ~ 2.3 star¬ 
forming regions are typically denser than local star¬ 
forming regions by an order of magnitude. We perform 
some tests to investigate the significance of the observed 
evolution in electron density. First, note that the SDSS 
line ratio distribution shown in Figure!^ is well-sampled 
and fairly narrow. We find that 89% of me SDSS sample 
has higher [S ii] ratios (lower electron densities) than the 
median z ^ 2.3 [S ii] ratio, while 64% of the z ^ 2.3 
sample has lower [S ii] ratios (higher electron densities) 
than the SDSS sample median. A two-sided Kolmogorov- 
Smirnov test on the z ^ 0 and z ~ 2.3 [S ii] distributions 


yields a K-S statistic of 0.452 and a p-value (the prob¬ 
ability that the two samples were drawn from the same 
underlying distribution) of 7.63 x 10“®, indicating that 
the z ~ 0 and z ~ 2.3 density distributions are signih- 
cantly different. 

We are attempting to estimate the electron density in 
HII regions, but [S il] emission coming from other com¬ 
ponents of the ISM could contaminate measured [S il] 
ratios from integrated-light galaxy spectra. The S ii ion 
has an ionization energy of 10.36 eV, somewhat lower 
than that of hydrogen. Because of its lower ionization en¬ 
ergy, the S II zone can extend beyond the boundary of an 
HII region. Additionally, [S ii]AA6716,6731 emission can 
be produced in diffuse ionized gas that is shock excited 
(Reynolds 1985 Martin|[l997 ). To investigate the effects 
of contamination from a dittuse ionized ISM component, 
we used measurements of [S ii] ratios of 44 HII regions 
in the star-forming spiral g alaxy NGC 628 observed as 
part of the CHAOS survey ( |Berg et ^|20I5 ). The spec¬ 
tra of these HII regions were attained by placing slits on 
top of individual HII regions and should contain very lit¬ 
tle light from the diffuse ISM. We find that the median 
[S ii]A67I6/A6731 ratio of these HII regions is 1.39, corre¬ 
sponding to an electron density of 38 cm“^. This ratio is 
nearly equivalent to the median [S ii] ratio of 1.41 for the 
SDSS sample, suggesting that our estimate of the typical 
local HII region density from SDSS is not significantly bi¬ 
ased by emission from diffuse ionized gas. However, all 
of these HII regions are from a single galaxy, and their 
median density may not be representative of the entire lo¬ 
cal HII region population. Currently, there is insufficient 
knowledge of the ISM structure of z ^ 2 galaxies to de¬ 
termine whether integrated-light spectra are signihcantly 
contaminated by emission from a diffuse component at 
that redshift. 

3.5. Electron density vs. galaxy properties 

We investigate whether the density of star-forming re¬ 
gions varies with other galaxy properties. As stated pre¬ 
viously, the nature of the function converting between 
line ratio and electron density makes it difficult to work 
with distributions in density space, especially when some 
objects have measured line ratios that are outside of the 
theoretically allowed values. For this reason we will look 
for relationships between density and galaxy properties 
using the line ratio as a proxy for the density. We plot the 
line ratios against stellar mass (M*), star-formation rate 
(SFR), and specific star-formation rate (sSFR; SF^M*) 
in Figure]^ The middle and right panels of Figures only 
include the subset of the z ^ 2.3 density sample with Ha 
and H/3 detections. We note that the z ~ 2.3 sample has 
signihcant overlap with the SDSS sample in M*, but the 
two are almost completely disjoint in SFR and sSFR. 
This difference is consistent w ith the evolution o f the 
SFR-M* relation with redshift ( Shivaei et al.|[2M5 ). We 
do not see evidence for any significant trends in line ra¬ 
tio (electron density) as a function of stellar mass, SFR, 
or sSFR among the local SDSS sample or the MOSDEF 
z ^ 2.3 sample. This observation is confirmed by per¬ 
forming a Spearman correlation test on each sample in 
each parameter space. No correlations are more signifi¬ 
cant than ~ 1(7. _ 

This result i s in co nflict with the recent work of IShi-l 
makawa et al. (2015), who found a correlation between 
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log(M,/MQ) log(SFR/M 0 yr-i) log(sSFR/Gyr-i) 


Figure 6. [S ll]AA6716,6731 (bottom row) and [O ll]AA3726,3729 (top row) as a function of stellar mass (left), SFR (middle), and sSFR 
(right). Blue and green points show the z ~ 2.3 [O ll] and [S ll] density samples, respectively. The gray two-dimensional histogram in the 
bottom row shows the distribution of the local comparison sample. Spectra from SDSS do not have a high enough spectral resolution to 
resolve the components of the [O ll] doublet. The blue, green, and gray dashed lines show the median line ratios for the corresponding 
sample. Dotted lines show the line ratios corresponding to denities of 10, 100, 1,000, and 10,000 cm“^. 


electron density and sSFR at the dcr level for 14 z ~ 2.5 
Ha emitters, and observed this correlation when stacking 
the spectra in two density bins. We do not see evidence 
for correlation among these properties with a larger and 
more repres entative sample that is generally consistent 


with that of 

Shimakawa et al. 

(2015 

) in stellar mass and 

SFR. 

Shima 

kawa et ai. ( 


found no correlation be- 


tween electron density and stellar mass, in agreement 
with our results. 


4. IONIZATION PARAMETER 

The ionization state of ionized gas in star-forming re¬ 
gions refers to the interaction between the ionizing source 
and ionized gas. This interaction modulates the relative 
populations of different ionic species, which directly influ¬ 
ence the observed emission line ratios. Multiple lines of 
evidence suggest that the ionization state in high-redshift 
star-forming regions may be systematically different from 
what is typically observed in the local universe in galax- 


ies of similar masses fe.g., iHolden et al. 112014 

Nakajima 

& Ouchi 

2014 

Steidel et aE 

2014|). i'he ionization pa- 


rameter is useful tor quantifying the ionization state, and 


encodes information about both the ionizing source and 
the surrounding ionized gas. In this section, we will use 
emission line measurements of z ^ 2.3 galaxies to inves¬ 
tigate the relationship between ionization parameter and 
other galaxy properties in order to probe the ionization 
state of gas in z ^ 2 star-forming regions. 

4.1. Definition of the ionization parameter 

We begin with some useful definitions related to the 
ionization state. Ionic species in HII regions are in ion¬ 
ization equlibrium, where the rate of ionization is equal 
to the rate of recombination. The ionization equilibrium 
condition for ionic species i can be written as 


where Ui is the number density of the ionic species, Qi is 
the rate of production of photons that can ionize species 
i, Ui is the effective ionization cross-section, is the 
number density of the once-ionized state of n^, rie is the 
electron density, and cxi is the recombination coefficient. 
If the ionization energy of Ui is close to that of hydro- 
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gen, which is the case for many of the ions that produce 
strong optical emission lines, then we can rearrange this 
expression and approximate the ratio of the relative pop¬ 
ulations in the higher and lower ionized states. 


nf ^ (jj Qo 

rii ai dTrr^rie 


( 10 ) 


Here, Qq is the rate of production of hydrogen-ionizing 
photons {hv > 13.6 eV). The first term (ai/ai) on the 
right hand side of this expression is a constant that 
will change for each ionic species. The second term 
(Qo/47rr^ne) contains only properties of the ionizing 
source and the gas, and is not dependent on the specific 
ionic species. The dimensionless ionization parameter U 
is defined as the second term on the right-hand side of 
equation [T0| divided by the speed of light, c. 


U = 


Qo 




( 11 ) 


Accordingly, the ratio of the relative population in an 
upper ionization state to that in a lower ionization state 
scales directly with the ionization parameter. Since the 
electron density is approximately the hydrogen gas den¬ 
sity in a fully-ionized plasma, the ionization parameter 
can be thought of as the ratio of the number density of 
hydrogen-ionizing photons to the number density of the 
hydrogen gas. When working with Hll regions it is con¬ 
venient to define the dimensionless ionization parameter 
using the radius of a canonical Stromgren sphere, Rs, as 
the distance between the gas and the ionizing source. 

U = - ( 12 ) 

‘i'KRgCne 

Often, the dimensional ionization parameter, q = c x U, 
is used instead, which is the ratio of the flux of ionizing 
photons at the Stromgren radius to the hydrogen number 
density. The definition of the Stromgren radius, based on 
a balance between ioniziation and recombination rates 
assuming case B recombination, is 

«, = ( , 13 ) 

\4:TraB£njj / \ATTaBen-^J 

where e is the volume filling factor of the gas. The volume 
filling factor can be defined by assuming that the gas 
is structured in dense clumps that are surrounded by a 
lower-density medium. In this case, the volume filling 
factor is defined as 


e = 


{ne^ 


(14) 


where (rie) is the global average electron density and rie.c 
is the electron density of the clumps. The volume fill¬ 
ing factor is equal to unity for a homogeneous constant- 
density gas, while its value decreases as the density of 
the clumps increases relative to the average density. We 
note that the density estimates from [S il] and [O ii] 
are based on luminosity-weighted measurements of emis¬ 
sion line strengths. Since emission strength scales as the 
square of the density, we are effectively measuring the 
clump density if a clu mpy gas geometry exists, not the 
global average density ( Kennicut^|1984 ). Using the def¬ 
inition of the Stromgreh sphere radius, we can simplify 


the ionization parameter and resolve its dependence on 
only the rate of ionizing photon production, the electron 
density, and the volume filling factor. 


U oc 


(15) 


The ionization parameter has a weak dependence on both 
the rate of ionizing photon production and the gas den¬ 
sity, and is somewhat more sensitive to the volume filling 
factor, through which this definition of the ionization pa¬ 
rameter contains information about the geometry of the 
gas. 

Defining the ionization parameter assuming the geom¬ 
etry of a Stromgren sphere is convenient, but likely does 
not hold for real HII regions. The Stromgren geometry 
assumes a sphere of constant density gas that immedi¬ 
ately surrounds the central ionizing source. In local HII 
regions, feedback from stellar winds can clear out a cav¬ 
ity around the ionizing star cluster such th at the ionized 
gas is a shell instead of a filled sphere (e.g., Watson et al. 
2008|). Accordingly, the Stromgren radius is not neces¬ 
sarily a good representative radius for the separation of 
the illuminated gas and the ionizing source. It is possi¬ 
ble that the wind-blown bubble geometry exists at high 
redshifts where the intensity of star formation is concen¬ 
trated, or some entirely different geometry such as inter¬ 
secting bubbles. The scalings presented in equation 15 
should then be used with caution because of the brea^I^ 
down of the Stromgren approximation. Real HII regions 
show a variety of complicat ed substructure and geometry 


(e.g., Pellegrini et al. 


20111. The Stromgren sphere defi- 


nition of the ionization parameter also assumes that the 
nebula is radiation-bounded (i.e. no hydr ogen-ionizing 


photons esca pe) instead of density-bounded (Nakajima & 


Ouchi 2014), which may not hold true at high redshifts. 


Additionally, in an integrated spectrum the measured 
ionization parameter is a luminosity-weighted average of 
all of the sources of emission inside the aperture, which 
includes multiple HII regions and emission from other 
ISM components. 

While the ionization parameter carries interesting in¬ 
formation about the ionizing source and gas geometry, 
it can be difficult to determine because it is not directly 
observable, but can only be estimated using calibrations 
derived from physically motivated models. The ioniza¬ 
tion parameter is often estimated using measurements 
of sets of emission line ratios that have some sensitiv¬ 
ity to the ionization state of the gas (e.g., lower- and 
higher-ionization states), in conjunction with the predic- 


20001 IKewiev & Dopita|20021 IDors et al. 112011 

Levesque 

Hi Richardson 

ditterences in i 

2014 

he t] 

Shirazi et al. 2014 
"anslation between o 

). Because of 

iservables and 


ionization parameter for different photoionization mod¬ 
els, it is convenient to instead use an empirical emis¬ 
sion line ratio as a proxy for the ionization parameter. 
Line ratios featuring both higher and lower ionization 
state transitions from the same element can be used to 
estimate the ionization parameter because of the rela¬ 
tion between the ionization parameter and the relative 
populations in the two ionization states. Here, we use 
032=[0 iii]AA5007,4959/[O ii]AA3726,3729 as a proxy 
for the ionization parameter. Systematic uncertainties 
in this approach result from the way in which a given 
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ionization parameter-sensitive line ratio depends also on 
the shape of the ionization spectrum and the metallicity 
of the gas, unless these two properties can be indepen¬ 
dently constrained. 

4.2. Sample selection 

In order to study the ionization state of high-redshift 
galaxies, we selected a sample of star-forming galaxies 
from the MOSDEF parent spectroscopic sample requir¬ 
ing objects to fall in the redshift range 2.0 < z < 2.6 
and have S/N> 3 in [O ii]AA3726,3729, H^, [O iii]A5007, 
and Ha. The flux of [O iii]A4959 is taken to be 1 /2.98 
of the [O iil]A5007 flux (|Storey & Zeippen||2000|). Ha 
and H/3 detections were necessary in order to correct line 
fluxes for dust attenuation, which is important for O 32 
because of the large wavelength separation of the emis¬ 
sion lines. Objects with Ha and H/3 detections also have 
robust dust-corrected SFRs based on Ha luminosities. 
AON were identihe d and removed ba sed on their X-ray 
and IR properties (Coil et al. 20151 and objects with 
[N ii]A6584 detected at ‘6a or greater were removed if 
log([N il]A6584/Ha) > —0.3. Any AON not removed by 
these selection criteria would introduce a bias in emis¬ 
sion line ratio diagrams. We discuss reasons why we are 
conhden t th at o ur s ample does not contain any AGN in 
Sections [T4| and [5^ Six additional objects were removed 
because of significant skyline contamination in the rele¬ 
vant emission lines. These criteria yield an ionization pa¬ 
rameter sample of 103 MOSDEF galaxies and the prop¬ 
erties of the sample are shown in Table This sample 
has properties that are nearly identical to those of the 
parent MOSDEF spectroscopic sample at z ~ 2.3, and is 
representative of star-forming galaxies with similar stel¬ 
lar masses at this redshift. 

In the discussion that follows, we will examine multi¬ 
ple line ratio diagrams, some of which involve additional 
emission lines along with oxygen and hydrogen Balmer- 
series strong lines. Therefore, we selected subsamples of 
the ionization parameter sample that have additional line 
detection criteria to plot in these spaces. We selected a 
subsample of 61 galaxies that additionally have S/N> 3 
in [N ii]A6585. We also selected a subset of 53 galaxies 
with S/N> 3 in both [N ii]A6585 and [S ii]AA6716,6731. 
The sample properties of the [N il] and [S ii] subsamples 
are presented in Table These stringent emission line 
cuts, requiring all or nearly all of the optical strong lines 
to be detected, introduce a bias against low-mass, low- 
SFR galaxies. However, the typical galaxy properties of 
the subsamples are not significantly different from those 
of the full ionization parameter sample or the MOSDEF 
parent z ~ 2.3 spectroscopic sample. While the [N ii] and 
[S ii] subsamples have slightly higher M* and SFR than 
the full sample, the median values are still consistent 
with the se galaxies falling on or near the z ^ 2 SF R-Ma, 
relation (Whitaker et al. |2014 Shivaei et al. 2015). 

We select a sample of typical star-forming galaxies 
in the local universe from SDSS DR7. We require 
0.04 < z < 0.1 and S/N>3 in [O ii]AA3726,3729, H/3, 
[O iii]A5007, Ha, and [N ii]A6585. Once again, the flux 
of [O iii]A4959 is assumed to be equal to 1/2.98 of the 
[O iii]A5007 flux dStorey & Zeippen||2000D. A detection 
in [N ii] is req uired because AGIN are r eiected using the 
demarcation of Kauffmann et al. (2003) in the [O iii]/H/3 
vs. [N ii]/Ha diagram. I'hese criteria yield a local sam¬ 


ple of 68,453 star-forming galaxies at z ~ 0.07. We also 
selected a subsample of 65,000 local galaxies that addi¬ 
tionally have S/N> 3 in [S ii]AA6716,6731. 

4.3. O 32 and global galaxy properties 

Many studies have suggested that z > 2 galaxies have 
systematically higher ionization parameters than are typ- 

ical for local galaxies (iHolden et al. 120141 iNakaiima et al. 

Kt 1 V i 'Ii joun n il ) I ' I lll Aun i l l . ■ 


20131 INakajima fc Uuchi||2014[ Smrazi et al. 2014[ btei- 
del et al.|2014p. However, it is imperative to consider the 
evolutioii in global galaxy properties with redshift when 
interpreting the apparently high ionization parameters 
observed at high redshifts. To this end, we investigate 
the dependence of the ionization parameter on global 
galaxy properties through the proxy of O. 32 , motivated by 


the cornparisons performed in Nakajima & Ouchi (2014). 
Figure presents O 32 vs. SFR and sSh'R for local star- 
forming galaxies (gray histogram) and the z ~ 2.3 ioniza¬ 
tion parameter sample (red circles), while Figurej^shows 
the dependence of O 32 on stell ar m ass for the same sam¬ 
ples. As described in Section [T5| the z ^ 2.3 galaxies 
span a similar range in M* as local SDSS galaxies, but 
have significantly higher SFR and sSFR at a given stel¬ 
lar mass, consistent with the evolution of the SFR-M^, 
relation with redshift ( Reddy et al.|2012 Whitaker et al. 
'2014 Shi^ei et al.||2015| . ' ' 

Figure]^ top panel, snows that there is a weak anti¬ 
correlation between O 32 and SFR for local galaxies, 
although there is signihcant scatter in O 32 at fixed 
SFR. Performing a Spearman correlation test yields a 
correlation coefficient of -0.19, indicating a weak anti¬ 
correlation, with a p-valuaj that is essentially zero due 
to the large sample size from SDSS. It is known that 
there is an anti-correlation between me tallicity and ion- 
ization parameter in the lo cal universe dPopita fc Evans 


19861 Dopitaet al.|2006a|b[|Perez-Montero|2U14[|SancIi^ 


et al. 2015), while SFR and metallicity a r e correlated 

with large scatter (Mannucci et al. 2010 Lara-Lopez 
|et al.||2013| ). Therefore, SFR and ionization parameter 
are anti-correlated (with large scatter) and the top panel 
of Figure [^confirms that O 32 acts as a proxy for the ion¬ 
ization parameter. The z ^ 2.3 galaxies also show a weak 
anti-correlation between O 32 and SFR, with a Spearman 
correlation coefficient of -0.30 and a p-value of 0.002. 

In the bottom panel of Figure!^ there is a tight cor¬ 
relation between sSFR and O 32 fOT z ^ 0 galaxies, with 
a Spearman correlation coefficient of 0.42. The bulk of 
local star-forming galaxies lie on this relation, while a 
small fraction of local galaxies with the least star for¬ 
mation (log(sSFR) < —1.0) do not show any correla¬ 
tion. It is possible that these low-sSFR galaxies have a 
very low level of ongoing star formation and are tran¬ 
sitioning to the red sequence via secular evolution or 
some other mechanism associated with the cessation of 
star formation. The ISM conditions in these low-sSFR 
galaxies may be different from those in typical local star¬ 
forming galaxies and would not be expected to follow 
the same trends. While z ~ 2.3 galaxies inhabit much 
of the same parameter space as highly star-forming local 

® In the Spearman correlation test, the p-value represents the 
probability of the dataset being drawn from an uncorrelated under¬ 
lying distribution. A correlation or anti-correlation with a p-value 
less than 0.003 has a significange greater than Scr. 
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Table 2 

Properties of the full ionization parameter sample, and [N ii] and [S ii] subsamples 



(z>- 



log(^) 

^^^0 med 

SFR** 

SFRmed^ 

Full sample 

2.29 

0.11 

8.97- 11.22 

10.0 

1.61 - 228 

23.8 

[N ii] subsample 

2.28 

0.11 

9.3 - 11.22 

10.22 

4.71 - 228 

33.2 

[S ii] subsample 

2.29 

0.11 

9.3 - 11.22 

10.26 

5.77- 228 

33.3 


^Average redshift of galaxies in the sample. 

^Standard deviation of the redshift distribution. 

^Range of log (M^/Mq) of galaxies in the sample. 

"^Median log(M*/MQ) of galaxies in the sample. 

®Range of SFR in M 0 /yr~^ of galaxies in the sample, determined from dust-corrected Hck luminosity. 
^Median SFR in M 0 /yr“^ of galaxies in the sample, determined from dust-corrected Ha luminosity. 
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Figure 7. O 32 vs. SFR (top) and sSFR (bottom) for local star¬ 
forming galaxies from SDSS (gray histogram) and z ~ 2.3 star¬ 
forming galaxies from MOSDEF (red circles). 




galaxies, they do not exhibit the same tight correlation 
between sSFR and O 32 . Performing a Spearman corre¬ 
lation test yields a correlation coefficient of 0.28 with a 
p-value of 0.004, indicating a weak but significant corre¬ 
lation. There is significant scatter in O 32 at fixed sSFR, 
and a non-negligible fraction of the high-redshift sample 


has high sSFR and low O 32 , a region of the parameter 
space where essentially no local galaxies are found. The 
difference between z ^ 0 and z ~ 2.3 galaxies in this 
space is interesting, but outside of the scope of this in¬ 
vestigation. 

The most intriguing of these diagrams is that of O 32 
vs. M*, shown in Figure!^ Local galaxies show a clear 
anti-correlation between O 32 and M* that is fairly tight, 
with a Spearman correlation coefficient of -0.52. In or¬ 
der to make the local trend more clear. Figure shows 
the running median O 32 at a given M* as a white line. 
This anti-correlation is consistent with the existence of 
a tight correlation between M* and me tallicity, known 
as the mass-metallicity relation (MZR; Lequeux et al. 

19791 [TVemonti etal. 120041 IKewley &: Fllison| 2 UU 8 | |Man- 

nucci et al.|2010 Andrews fc Martini|2013p, and an anti- 


correlation between metallicity and ionizat ion parameter 
( Perez-Montero|[2014 Sanchez et al.|[20l5 ). In fact, O32 


has been show n to be a metallicity indicator for object s 
up to z ^ 0.8 (|Maiolino et al.||2008t | Jones et al.||2015|). 


The relation between U 32 and Ivi* appears to flatten out 
at high stellar masses, consistent with the observed be¬ 
havior of the local MZR. 

The z ~ 2.3 sample also displays a fairly tight relation 
between O32 and M*such that higher M* corresponds 
to lower O32, with a Spearman correlation coefficient of 
-0.57 and a p-value of 3.7 x 10“^°. We show the run¬ 
ning median O 32 of the high-redshift galaxies at a given 
M* as a dark red line in Figure The high-redshift 
anti-correlation shows nearly the same slope as that of 
the local relation, only offset towards higher O32 at fixed 
M*. We find that z ^ 2.3 galaxies have O32 values ~ 0.6 
dex higher at a given M*, suggesting that high-redshift 
galaxies have significantly higher ionization parameters 
than local galaxies of the same stellar mass if the trans¬ 
lation between O32 and ionization parameter is the same 


lated at z ~ 2.3 as well fe.g., 

Erb et al. 2006 

Maiolino 

et al. 

2008 

Steidel et al.||2U14 Sanders et al.||2 

OTS), but 


the MZR evolves such that galaxies at a given iVR 
lower metallicities than are observed locally. The exis¬ 
tence of a clear O32 vs. M* anti-correlation at z ^ 2.3 is 
suggestive of an anti-correlation between metallicity and 
ionization parameter existing at high redshifts as well. 

The striking similarity of the shape of the O32 vs. M* 
relation for local and z ~ 2.3 galaxies suggests that a 
similar mechanism may set the observed ionization pa¬ 
rameter at both redshifts, but must evolve with redshift 
such that high-redshift galaxies have higher ionization 
parameters at a given stellar mass. The evolution of the 
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Figure 8 . O 32 vs. stellar mass for local star-forming galaxies from SDSS (gray histogram) and z ~ 2.3 star-forming galaxies from 

MOSDEF (red circles). The solid white line shows the running median O 32 in bins of stellar mass for local galaxies. The solid dark red 
line shows the running median of z ~ 2.3 galaxies. 


MZR from z ^ 0 to z ~ 2.3 provides a natural expla¬ 
nation for the apparent change in ionization parameter 
at fixed stellar mass. We further investigate the inter¬ 
play of O32, M*, and metallicity by employing the use of 
emission line ratio diagrams. 

4.4. O32 and metallicity 

We have shown that there is a relation between O32 and 
M*, the shape of which is similar for local and z ~ 2.3 
galaxies, that has a higher normalization in O32 at high 
redshift. However, stellar mass is a global property that 
is not directly related to the production of emission lines 
in individual star-forming regions that are an observable 
probe of the ionization parameter. The metallicity of 
the gas in star-forming regions, on the other hand, has a 
direct impact on both the ionizing spectrum, assuming 
the gas-phase metallicity is related to the stellar metal¬ 
licity of the ionizing cluster, and the intrinsic emission 
line fluxes. Therefore, comparing O32 values at fixed 
metallicity rather than fixed M* utilizes a property that 
directly influences the physical conditions in star-forming 
regions, including ionization parameter, and removes sys¬ 
tematic effects introduced by the evolution of galaxy scal¬ 
ing relations with M*. 

Figurej^shows the dependence of O32 on R23 (left) and 
03N2 (ri^t) for local SDSS galaxies (gray histogram) 
and the z ~ 2.3 ionization parameter sample. Note that 
the right-hand panel only includes high-redshift galaxies 
in the [N ii] subsample since the 03N2 ratio requires 
[N ii]A6584. Figure 0 demonstrates that the z ~ 2.3 
sample is likely free from AGN contamination. AGN 
are found at low 03N2, high R23, and high O32. No 
galaxies in the z ~ 2.3 sample fall on the AGN sequence 


in the 03N2 diagram. In the R23 diagram, there are 
several z ~ 2.3 galaxies in the same region of parameter 
space as the AGN sequ ence (log(R23) > 0 . 95), b ut these 
galaxies fall below the Kauffmann et al. (2003) line in 
the [O iii]/H/3 vs. [N Iij/Ha diagram and tend to have 
large uncertainties in R23. 

The R23 index is sensitive to m etallicity, but is double¬ 
valued (Kewley & Dopita 2002). This problem can be 
overcomie by using a secon d excitation-sens i tive li ne ra¬ 
tio in tandem with R23. |Shapley et al. (2015) used 


direct-method abundan ces of stacks of local galaxies (An¬ 


drews & Martini 2013) to show that the local sequence 


in the (J 32 vs. R 23 diagram is a sequence in monotoni- 
cally increasing metallicity from the high-032, high-R23, 
high-excitation tail towards the I 0 W-O 32 , I 0 W-R 23 , low- 
excitation region. The 03N2 index is also a metallicity 
indicator, reflecting the fact that the narrow sequence 
of star-forming galaxies in the BPT diagram is also a 
sequence in metallicity, and has been used as such in 


empir ical metallicity calibrations (e.g., Pettini & Pagel 


20041. The position on either of these diagrams reflects 


the ionization parameter at a given oxygen abundance. 

Using a smaller sample f rom the early MOSDEF 
dataset, |Shapley et al.| (|2015D showed that z ~ 2.3 star¬ 
forming galaxies seem to follow the same distribution as 
local galaxies in the low-metallicity, high-excitation tail 
of the O32 vs. R23 diagram. We confirm this finding 
with the current MOSDEF z ~ 2.3 sample. The z ~ 2.3 
galaxies display no systematic offset with respect to the 
local galaxies, and the bulk of the high-redshift sample 
inhabits the region in which 12 -|- log(0/H) < 8.6 for 
local galaxies. The running median of the z ~ 2.3 sam¬ 
ple (dark red line) closely follows that of z ~ 0 galaxies 
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Figure 9. O 32 vs. R 23 (left panel) and 03N2 (right panel). The gray histogram shows the distribution of local star-forming galaxies from 
SDSS. The red points and error bars denote 2 ~ 2.3 star-forming galaxies from the MOSDEF survey. In both panels, the white line shows 

the running median of the local sample. The dark red line shows the running median of the z ~ 2.3 in the left panel, and the best-fit linear 

relation to the 2 ~ 2.3 galaxies in the right panel. The right panel only includes galaxies in the [N ii] subsample. 

(white line) in the O32 vs. R23 diagram]^ We also find 
that z ~ 2.3 galaxies closely follow the subsolar abun¬ 
dance tale of the local distribution in the O32 vs. 03N2 
diagram, although there is a slight systematic offset to¬ 
wards lower 03N2 at fixed O32. By comparing the best- 
fit linear relation of the z ~ 2.3 sample (dark red line) to 
the running median of the z ^ 0 sample (white line), we 
find the mean value of this offset to be 0.23 dex in 03N2 
at fixed O32. We will di scuss the offset in the 03N2 
diagram further in Section [5?^ 

We propose the following scenario to explain the ob¬ 
served position of high-redshift galaxies in the O32 vs. 

R,23, O32 vs. 03N2, and O32 vs. M* diagrams. Taken 
together, the two plots in Figure [^suggest that the ion¬ 
ization parameter at fixed metalTicity is the same for 
galaxies at z ^ 2.3 and z ~ 0. Galaxies at high red- 
shift must then follow the same anti-correlation between 
oxygen abundance and ionization parameter that is ob¬ 
served locally. This U vs. 0/H relation shows little to no 
redshift evolution between z ^ 0 and z ~ 2.3. The nearly 
constant offset towards higher O32 at fixed M* observed 
in Figure^ is then simply a co nsequence of the evolu tion 
of the M2 
found z ^ 


with redshift. In Sanders et al. (2015), we 
2.3 galaxies have metallicities ~ 0.3 dex 


ower 


than lo cal galaxies at fixed M» using the 03N2 calibra¬ 
tion of|Pettini & Pagel| (|2004|. If z ~ 0 and z ~ 2.3 
galaxies follow the same metallicity-ionization parame¬ 
ter relation, then the decrease in metallicity at fixed M* 
leads to an increase in O32, as observed in Figure]^ 
Earlier results suggesting that high-redshift g^xies 
have higher ionization parameters than are seen locally 
were based on either a c omparison at fixed stellar mass 
(e.g., Holden et al.||2014 | or a comparison of the average 
ionization parameter of the entire local star-forming pop- 


® A running median creates a smooth representation of the local 
sequence of star-forming galaxies because of the large number of 
galaxies in the local sample. The small number of galaxies in the 
2 ~ 2.3 sample and subsamples can cause the running median 
to be uneven and erratic. A linear fit better represents the 2 ~ 
2.3 subsample distributions in the 03N2 diagram, while a running 
median is still used in the R 23 diagram because the O 32 vs. R 23 
sequence displays curvature. 


ulation to that of high-redshift galaxies (e.g., Hainlinej 
et al.||2009). At fixed metallicity, which is more directly 
related to the ionization state of the gas, we find that 
z ~ 2.3 galaxies have roughly the same ionization param¬ 
eters compared to local galaxies. The ionization state of 
star-forming regions in galaxies in our high-redshift sam¬ 
ple must then be similar to what is observed in metal- 
poor local galaxies. 

In Section we found that the density of star-forming 
regions increases significantly from z ^ 0 to z ^ 2.3. 
Here, we are suggesting that galaxies at z ~ 0 and 
z ~ 2.3 have the same ionization parameter at fixed 
metallicity despite the difference in density. This ini¬ 
tially seems to be at odds with the scaling of ionization 
parameter with density presented in equation |15[ An 
increase in density of an order of magnitude would cor¬ 
respond to an increase in the ionization para meter of 
more than a factor of two. Shirazi et al. (2014| used the 
relation between ionization parameter and electron den¬ 
sity to explain the high ionization parameters observed 
at z ~ 2 — 3 by an increase in the density of star-forming 
regions. If the scalings in equation [T5] hold at z ^ 0 and 
z ~ 2.3, then an increase in density does not guarantee 
an increase in the ionization parameter unless the volume 
filling factor is the same at both redshifts. If high-redshift 
HII regions are clumpier, the volume filling factor would 
decrease with redshift and could offset the effect of an 
increase in density. However, the scalings presented in 
equation]^ are derived assuming the Stromgren approx¬ 
imation, wTiich does not apply to many local HII regions 
and likely does not hold at z ~ 2.3. Thus, an increase of 
a factor of 10 in density from z ~ 0 to z ~ 2.3 does not 
necessitate higher ionization parameters at fixed metal¬ 
licity. 

Our proposed scenario is only valid if some assump¬ 
tions hold true. A location in O32 vs. R23 space must 
correspond to the same metallicity regardless of redshift. 
The relation between 03N2 and metallicity must not 
evolve significantly with redshift. Finally, the translation 
between ionization parameter and O32 must not evolve 
significantly from z ~ 0 to z ^ 2.3. We investigate the 
validity of these assumptions in the following section. 
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5. DISCUSSION 

In this section, we present evidence in support of our 
proposed scenario that z ^ 2.3 star-forming galaxies have 
roughly the same ionization parameter as z ^ 0 galax¬ 
ies when comparing at fixed metallicity. In Section |5.1 
we investigate whether the translation between O 32 anc 
ionization parameter evolves with redshift using a set of 
simple photoionization models in combination with the 
observed position of z ^ 2.3 galaxies in emission line ra¬ 
tio diagnostic diagrams. We explore whether or not there 
is evidence that the metallicity dependence of R 23 and 
03 N2 ch anges signif icant ly from z ~ 0 to z ^ 2.3 in Sec¬ 
tion 5.2[ In Section |5.31 we discuss additional evidence 
for the existence of an anti-correlation between O/H and 
U and how the existence of this anti-correlation affects 
the interpretation of photoionization model grids. In Sec¬ 
tion |5^ we use our results to explore the cause of the 
well-known offset of z ~ 2 galaxies in the [N ii] BPT dia¬ 
gram compared to local galaxies. Finally, we discuss the 
uncertainty that diffuse ionized gas introduces in the in¬ 
terpretation of emission line ratios from integrated-light 
galaxy spectra in Section [57^ 


5.1. Is the ionization parameter-O 32 relation redshift 

invariant? 

We first address whether or not the relationship be¬ 
tween ionization parameter and O 32 is the same at z ~ 0 
and z ^ 2.3. We describe the simple photoionization 
models used for this analysis, discuss the interdepen¬ 
dence of the shape of the ionizing spectrum, ionization 
parameter, and O 32 , and estimate the typical ionization 
parameter of z ~ 2.3 and local galaxies based on the 
models. 

5.1.1. Description of the Cloudy photoionization models 

As previously mentioned, the chief difficulty in de- 
terming the ionization parameter is that it can only be 
done with reference to a specific set of photoionization 
models. The extent to which the value of the estimated 
ionization parameter can be trusted depends on how 
well the models represent the observed objects and can 
produce self-consistent predictions in multiple line ratio 
spaces. Given the uncertainty of photoionization mod¬ 
els both locally and at high-redshift where the physical 
properties are less constrained, we use a suite of simple 
photoionization models to understand qualitatively the 
trends in emission line ratios when changing the differ¬ 
ent input parameters of the models. We do not, however, 
use these models to place tight constraints on the metal¬ 
licity or ionization parameter of any local or high-redshift 
galaxies. 

We use the photoionization code Cloudj|^ to model 
emission line ratios from star-forming regions with a 
range of physical con ditions. These model s are very sim¬ 
ilar to those used by Steidel et al. (2014) to investigate 
the position of z ~ 2.3 galaxies in the [O iil]/H/3 vs. 
[N ii]/Ha diagram. There are five main input parame¬ 
ters that determine the location of a grid point in vari¬ 
ous emission line ratio diagrams: hydrogen gas density, 
gas-phase metallicity, ionization parameter, shape of the 
ionizing spectrum, and N/0 abundance ratio. 

Calculation s were performed wi th version 13.02 of Cloudy, 
last described by|Ferland et aL|(|2013[|. 


In Section!^ we presented a robust characterization of 
the typical electron density in local and z ~ 2.3 galaxies, 
finding densities of ~ 25 cm“^ and ^ 250 cm“^, respec¬ 
tively. In HII regions, the gas is fully ionized and the elec¬ 
tron density provides a good estimate of the hydrogen gas 
density. Since we have measured the typical density at 
z ^ 0 and z ~ 2.3, we only allow the density to be either 
25 cm“^ or 250 cm“^ in the models. The metallicity sets 
the global abundance scale of the gas, which is assumed 
to follow a solar abundance pattern with the exception 
of nitrogen. We vary the metallicity between 0.2 Zq and 
I.O Zq (12-|-log(0/H) = 8.0 — 8.69) in 0.2 Zq steps with 
solar metallicity corr esponding to 12 -|-log(O/H )0 = 8.69 
(Asplund et al.|2009). The ionization parameter sets the 


ionization state ot the gas and is allowed to vary between 
\og{U) = —3.6 and —1.5 (log(^;j^) = 6.9 — 9.0) in 0.1 

dex steps. _ 

As described in Section |4.1[ the shape of the ionizing 
spectrum affects the relative populations of an element in 
different ionized states. A harder ionizing spectrum re¬ 
sults in a larger fraction of oxygen in O ill compared 
to O II, for example. We use a blackbody spectrum 
with an effective temperature of 40,000 K, 50,000 K, or 
60,000 K as the input ionizing spectrum. Steidel et al. 
(2014) showed that, when properly normalized, a black- 
body spectrum is a good approximation of the spectrum 
of massive stars bluewards of 912 A using BPASS stellar 
models that include effects from b inarity (Eldridge et al. 


2011 Eldridge & Stanway 2012). We have found this 
observation to hold true when using Starburst99 (SB99) 
stellar models that include effects of rotation in massive 


stars (Leitherer et al. 2014). We note that the effective 
temperature ot the blackbody is not the same as the ef¬ 
fective temperature of a star. It is simply a parameter 
that allows us to specify the shape of the input spec¬ 
trum. When referring to a “harder” ionizing spectrum, 
we are referring to an increase in the blackbody effective 
temperature of the input spectrum. We additionally uti¬ 
lize input spect ra produced by SB99 using the Geneva 
2012/13 tracks (Ekstrom et al. j2012 Georgy et al.|2 0l'^ 
that include the effects ot rotation in massive stars dLei- 


therer et al. 2014). We create two input spectra from 
bBOO that create bracketing cases of a very hard ioniz¬ 
ing spectrum, which we refer to as ^^SB99 hard" , and 
a softer ionizing spectrum, which we refer to as ^^SB99 
soft.” SB99 hard is produced assuming a single burst 
of star-formation that formed 0.5 Myr ago with stel¬ 
lar metallicity of 1/7 Zq. SB99 soft instead assumes 
a 10 Myr-old population with solar metallicity formed 
with a continuous SFR of 1 MQ/yr. An age of 10 Myr 
was chosen to ensure that the ionizing spectrum of the 
stellar populatio n had reached a ste ady state, occuring 
after ~ 5 Myr (Kewley et al. 2001). In both cases, a 


Salpeter (1955) liVlh' slope is assumed above 0.5 M©. 
'I'hese two cases roughly bracket the range of ionizing 
spectra appropriate for HII regions contributing signifi¬ 
cantly to integrated-light galaxy spectra. 

In reality, the shape of the ionizing spectrum should 
be related to the gas-phase metallicity, which traces the 
stellar metallicity since recently formed massive stars 
are ionizing the remnants of their birth cloud. In lower 
metallicity stars, there is less metal line blanketing and 
opacity in the stellar atmospheres, leading to hotter ef- 
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fective temperatures and harder ionizing spectra. How¬ 
ever, we allow the shape of the ionizing spectrum to vary 
separately from the metallicity to accomodate the possi¬ 
bility that the hardness of the ionizing spectrum at fixed 
metallicity evolves with redshift. 

A solar abundance pattern is assumed for all elements 


izing spectrum evolving with redshift by combining the 
observed line ratios for our z ~ 2.3 sample with grids 
from the Cloudy models in emission line ratio diagrams. 
It has been proposed that the position of high-redshift 
galaxies in the [O iii]/H/3 vs. [N ii]/Ha diagram can be 
explained by a systematically harder ionizing spectrum 


dance ratio is observed to have a dependence on 0/H, licity ( 

Kewley et al. 

|2013a 

Steidel et al. 

2014 

). If true, 

such that IN/(J is a constant value at low abundance but this explanation would leac 

to evolution ot the ionization 


begins to rise roughly linearly with 0/H at higher abun¬ 
dance 2009J_ 


2012 An5rews°^°^arnm ' 2L)13[ ~ FereS^MoSCTO ' 2U14p . 


At low metallicity, nitrogen is a primary nucleosyn- 
thetic product of hydrogen and helium burning. At high 
metallicity, nitrogen is produced as a secondary product 
through the CNO cycle where the yield of nitogren de¬ 
pends on the amount of pre-existing C and O, which leads 


to th e dependence of N/0 on 0/H (|van Zee & Haynes 
2006 1 . There is disagreement ab out the shape of the IN /(J 


vs. (J/H relation (see Fig. 12 in iSteide l et al.|2014 |. We 
assum e the relation found by |Herez-Montero fe Contini 


( 20091) , which is a simple linear relation over the range of 


metaliicities considered in the models. We note that the 
assumed N/0 ratio can strongly affect those line ratios 
involving [N ii]A6584, but has negligible effects on other 
line ratios. 

5.1.2. Ionization parameter and the hardness of the ionizing 
spectrum at z ^ 0 to z ~ 2.3 

From the set of models with density, gas-phase metal¬ 
licity, ionization parameter, ionizing spectrum, and N/0 
ratio defined as above, we extract the relationship be¬ 
tween O 32 and ionization parameter to resolve its de¬ 
pendence on the various input parameters. Figure |10| 
shows O 32 vs. log(U) for models with electron density m 
25 cm“^ (top panel) and 250 cm“^ (bottom panel). The 
line color denotes the input spectrum, while the different 
lines of a single color connect grid points with the same 
metallicity. We note that the SB99 hard models behave 
very similarly to the 60,000 K blackbody models. The 
SB99 soft spectrum appears to be slightly softer than a 
50,000 K blackbody based on the position of the model 
grids. Variation in N/0 has no effect on the relationship 
between O 32 and U. In addition. Figure IT shows that 
this relation has very little dependence onthe gas den¬ 
sity. Lines of constant metallicity show that the relation 
between O 32 and U has only a small dependence on gas- 
phase metallicity when the shape of the input spectrum 
is fixed. The one exception is for a soft ionizing spec¬ 
trum with high ionization parameter and low gas den¬ 
sity (Teff = 40,000 K, \og{U) > —2.5, Ue = 25 cm“^), a 
region in which real objects are unlikely to be found. Cal¬ 
ibrations of O 32 and the ionization parameter typically 
show significant dependence on metallicity because the 


ionizing spectrum is tied to the metallicity (e.g., Kew-, 
ley fc Dopita||2002 ). On the other hand, the hardness of 


the ionizing spectrum has a significant effect on the O 32 
vs. logilA) relation, such that a harder ionizing spectrum 
produces a larger O 32 value at fixed ionization parame¬ 
ter. Therefore, the translation between O 32 and ioniza¬ 
tion parameter will only show significant evolution if the 
hardness of the ionizing spectrum at a given metallicity 
evolves with redshift. 

We investigate the possibility of the shape of the ion¬ 


paramet er -0 32 relation. 

Figure 11 shows the [O iii]/H/3 vs. [N il]/Ha ([N il] 
BPT) diagram (left column) and the [O iil]/H,d vs. 
[S ii]/Ha ([S ii] BPT) diagram (right column). The lo¬ 
cal distribution of star-forming galaxies is shown as the 
gray histogram in all panels. The top row presents the 
observed line ratios for the z ~ 2.3 [N ii] subsample (top 
left) and [S ii] subsample (top right). The middle and 
bottom rows show the Cloudy model grid points as cir¬ 
cles. The size of the circle indicates the metallicity, with 
the largest size indicating solar metallicity and the small¬ 
est size indicating 0.2 Z©. The color indicates the input 
ionizing spectrum, with the effective temperature listed 
for blackbodies in the middle row and the bracketing 
SB99 models shown in the bottom row. Solid lines con¬ 
nect points of constant ionization parameter, with ioniza¬ 
tion parameter increasing to the upper left. The model 
grids are only shown for a density of 250 cm“^. Dis¬ 
playing models at a single density will suffice since we 
are only interested in discerning trends in line ratio with 
the model input parameters instead of making quantita¬ 
tive predictions. While the absolute line ratios change, 
the trends with metallicity, ionization parameter, and 
hardness of the ionizing spectrum are the same regard¬ 
less of the assumed density. Qualitatively, increasing the 
electron density while keeping all other parameters fixed 
increases [O iil]/H/3, [N il]/Ha, and JS iij/Ha. 

In the top right panel of Figure [l^ we observe the 
well-documented offset of z ^ 2.3 galaxies towards higher 
[O iii]/H/3 and/or [N ii]/Ha compared to the local star¬ 
forming sequence. The magnitude of this offset can be 
observed by comparing the running median of the z ~ 2.3 
sample (dark red line) to that of the z ^ 0 sample (white 
line). However, we do not observe a significant offset 
between z ~ 0 and z ~ 2.3 galaxies in the [S ii] BPT 
diagram (top left pan el), in agreement w it h results from 
early M OSDEF data (Shapley et al. 2015). Masters et al. 
(2014) found a similar result using a composite spectrum 
ot 24 z ^ 2 emission-line galaxies from the WISP survey. 
In the [S ii] BPT diagram, the running median of the 
z ~ 2.3 sample is marginally offset to lower [O iii]/H/3 
and/or [S Iij/Ha, although this offset is not significant 
given the sample size and measurement uncertainty. In 
the middle and lower left panels, it can be seen that an in¬ 
crease in the hardness of the ionizing spectrum generally 
moves grid points to higher [O iii]/H/3 and [N ii]/Ha at 
fixed metallicity and ionization parameter, which could 
potentially explain the z ~ 2.3 offset. In the [S ii] BPT 
diagram (middle and lower right), the models show that 
an increase in the hardness of the ionizing spectrum in¬ 
creases [O iii]/H/3 at fixed [S Iij/Ha. We emphasize that 
we are using these models to demonstrate how predicted 
line ratios change qualitatively as the hardness of the 
ionizing spectrum varies. We are less concerned with 
a match between a specific blackbody or SB99 spectrum 
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Figure 10. O 32 vs. ionization parameter, ti, from simple photoionization models of HII regions assuming a gas density of 25 cm ® (top 
panel) and 250 cm“® (bottom panel). The dependence of O 32 on lA is extracted from a set of models assuming different ionizing spectra, 
denoted by the color of the line, and different gas-phase metallicities. Input ionizing spectra are assumed to be either a blackbody spectrum 
with an effective temperature of 40,000-60,000 K or one of two spectra produced by Starburst99 {SB99 hard and SB99 soft). Each line of 
a single color connects models with the same gas-phase metallicity, from 0.2 Zq to 1.0 Zq. 


and either the local or high-redshift excitation sequences. 

In both BPT diagrams, an increase in the hardness of 
the ionizing spectrum leads to a signihcant increase in 
[O iil]/H/3 at fixed [N Ilj/Hct or [S Iij/Ha for the same 
ionization parameter and metallicity. If a harder ioniz¬ 
ing spectrum was the cause of the offset of high-redshift 
galaxies in the [N ii] BPT diagram, simple photoioniza¬ 
tion models predict that there should also be an offset 
of similar magnitude towards higher [O Iii]/H,0 in the 
[S ii] BPT diagram. We observe no significant offset be¬ 
tween the z ^ 2.3 and z ^ 0 galaxies in the [S ii] BPT 
diagram. We note that in the [N il] BPT diagram, the 
60,000 K and 50,000 K blackbody models approximately 
match the position of the z ^ 2.3 and local galaxies, re¬ 
spectively. The same 60,000 K and 50,000 K blackbody 
models are offset from one another in the [S ii] BPT dia¬ 
gram, and yet no offset is observed in the corresponding 
observations of z ^ 2.3 and local galaxies. We conclude 
that there is not a systematically harder ionizing spec¬ 
trum in high-redshift star-forming regions compared to 
local star-forming regions of similar metallicity. We will 
revi sit th e cause of the [N ii] BPT diagram offset in Sec¬ 
tion |5.4| Having established that the relation between 
ionization parameter and O 32 is only significantly sensi¬ 
tive to changes in the hardness of the ionizing spectrum, 
we further conclude that the U-O 32 relation does not 
strongly evolve with redshift. 

5.1.3. Estimating the typical ionization parameter at z ^ 2.3 

Having shown that the translation between O 32 and 
ionization parameter does not appear to evolve signifi¬ 
cantly with redshift, we will now estimate the range of 
ionization parameters that would be inferred from the lo¬ 
cal and z ~ 2.3 galaxy samples. We advise caution when 


using these ionization parameter estimates because a dif¬ 
ferent ionizing spectrum should be used for the low- and 
high-metallicity galaxy populations and the use of differ¬ 
ent stellar models will change these estimates. Matching 
the appropriate stellar models to a given object is fur¬ 
ther complicated by uncertainty in the metallicity esti¬ 
mates used to tie the observed metallicity to the absolute 
metallicity built into the stellar models. Since we do not 
know which input ionizing spectrum is appropriate for 
each sample, we provide estimates of the ionization pa¬ 
rameter assuming each of the five ionizing spectra con¬ 
sidered in the models. The median log( 032 ) value of 
the z ^ 2.3 sample is 0.10 and the middle 68 % span 
log( 032 ) = —0.11 to 0.37. The local SDSS sample has 
a median log( 032 ) value of —0.53 and the middle 68 % 
span log( 032 ) = —0.71 to —0.24. Estimates of the corre¬ 
sponding values of log(^) are presented in Table [fusing 
the curves in Figure [l^ We note that the median ioniza¬ 
tion parameter is approximately the same for z ^ 0 and 
z ^ 2.3 galaxies if the ionizing spectrum of local galaxies 
is well-described by a 40,000 K blackbody and that of 
z ~ 2.3 galaxies is described by a 60,000 K blackbody. 
However, we have shown that z ~ 0 and z ~ 2.3 galaxies 
have similar ionizing spectra at fixed metallicity. While 
the ionizing spectrum at the median metallicity of each 
sample will be different, it is unlikely that the magnitude 
of that difference is as large as the difference between a 
40,000 K and 60,000 K blackbody spectrum. Therefore, 
z ~ 2.3 galaxies have a higher median ionization param¬ 
eter than local galaxies. 


5.2. Is the dependence of R 23 and 03N2 on metallicity 
redshift invariant? 
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Figure 11. [N ii] (left column) and [S ii] (right column) BPT diagrams. The local sequence of star-forming galaxies from SDSS is shown 
as the gray histogram in all panels. The top row shows the position of 2 ~ 2.3 star-forming galaxies in the [N ii] and [S ii] BPT diagrams 
as red points with error bars. The white and dark red lines show the running median of the local and 2 ~ 2.3 samples, respectively. The 
bottom two rows show grids from simple photoionization models of HII regions produced using Cloudy. Different colors represent different 
input ionizing spectra. In the middle row, the input ionizing spectrum is assumed to be a blackbody spectrum and the effective temperature 
is listed, representing the hardness of the ionizing spectrum. In the bottom row, the input ionizing spectrum is assumed to be one of two 
spectra produced by the stellar population synthesis code Starburst99 {SB99 hard and SB99 soft). Colored circles show the model grid 
points, where the size of the circle represents the gas-phase metallicity, from 0.2 Zq( smallest circles) to solar metallicity (largest circles). 
Grid points of constant ionization parameter with the same input ionizing spectrum are connected with solid lines. All model grids shown 
are calculated assuming a hydrogen gas density of 250 cm~^. 
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Table 3 

Ionization parameter estimates of local SDSS and .2 2.3 star-forming galaxies based on O32 for each of the five input ionizing spectra 

assumed in our photoionization models. 



log(032) 

log(7y) assuming an ionizing spectrum of: 

40,000 K 50,000 K 60,000 K SB99 soft SB99 hard 


median 

-0.53 

-2.90 

-3.20 

-3.40 

-3.20 

-3.45 

SDSS 

lower 68%“^ 

-0.71 

-3.05 

-3.35 

-3.50 

-3.35 

-3.60 


upper 68% 

-0.24 

-2.60 

-3.05 

-3.20 

-3.00 

-3.25 


median 

0.10 

-2.10 

-2.75 

-2.95 

-2.70 

-3.00 

2 2.3 

lower 68%“^ 

-0.11 

-2.45 

-2.95 

-3.10 

-2.90 

-3.15 


upper 68%’’ 

0.37 

^ -1-5'' 

-2.50 

-2.70 

-2.40 

-2.75 


^The lower bound on the middle 68% of the distribution of log(032) 

'^The upper bound on the middle 68% of the distribution of log(032) 

^The range of ionization parameters considered in this set of models did not extend high enough to give log(032) = 0.37 for a 40,000 K 
blackbody input spectrum, so a lower limit is assigned. 
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We proposed that the relationship between metallic- 
ity and ionization p ara met er is the same locally and at 
z ~ 2.3. In Figures [T0| and 0 we have shown evidence 
that the relationship between ionization parameter and 
O 32 does not significantly change with redshift. In or¬ 
der for the proposed scenario to hold in concordance 
with Figure]^ we must also show that the dependence 
of R 23 and 03N2 on metallicity does not significantly 
evolve with redshift. In Figure we show the local 
star-forming sample (gray histogram) and the z ^ 2.3 
[N ii] subsample (red circles) in the spaces of O 32 vs. 
R 23 (left), 03N2 (middle), and N2 (right). 

As previously mentioned, no systematic offset is ob¬ 
served between z ~ 0 and z ^ 2.3 galaxies in the O 32 
vs. R 23 diagram, although z ^ 2.3 objects only occupy 
the low-metallicity tail of the local distribution. We also 
pointed out the 0.23 dex offset z ~ 2.3 g alaxi es show to¬ 
wards low 03N2 at fixed O 32 in Section [F4l We find a 
larger systematic offset between local and high-redshift 
galaxies in the O 32 vs. N 2 diagram of 0.33 dex higher 
N2 at fixed O 32 . There is an assymetric scatter towards 
high N2, such that no z ~ 2.3 galaxies scatter below the 
local sequence. We additionally observed no significant 
offset in the [S ii] BPT diagram (Fig. II). These results 
collectively suggest that, on average, hipi-redshift galax¬ 
ies have higher [N ii]A6584 flux compared to the strength 
of other strong optical emission lines than typical local 
galaxies with the same ionization parameter. Outside of 
line ratios involving nitrogen, high-redshift galaxies ap¬ 
pear to behave similarly to the low-metallicity tail of the 
local distribution. We have measurements constraining 
the evolution of density with redshift, and have found no 
evidence suggesting the ionization parameter or hard¬ 
ness of the ionizing spectrum evolve significantly at fixed 
metallicity. Therefore, we conclude that metallicity in¬ 
dicators involving nitrogen will evolve with redshift, but 
the relation between metallicity and location in the O 32 
vs. R 23 diagram is likely the same at z ~ 0 and z ~ 2.3. 
Jones et al. (2015) recently found that the relationship 
between R 23 , U 32 , and di rect-method met allici ty does 


not evolve out to z ^ 0.8. Liu et al. (2008) and 


et al. (2014) concluded that the IN 2 indicator significantly 


Steidel 


over-estimates the metallicity of objects offset from the 
local star-forming sequence, while the 03N2 indicator 
has a much smaller bias. In co njunction with the argu¬ 
ments presented in Section [O] these results support our 
proposed scenario that galaxies at fixed metallicity have 
the same ionization parameter locally and at z ~ 2.3. 

5.3. Additional evidence for the existence of the 
metallicity-ionization parameter anti-correlation 

The argument for an anti-correlation between 0/H 
and ionization parameter arises from the fact that low- 
metallicity stars produce more ionizing photons in to¬ 
tal an d have harder ionizing spectrum (Leitherer et al. 
2014). Observational evidence for this anti-correlation 
comes from finding tha t O 32 and [O ill]/H/9 increase with 
decre asing metallicity (Maiolino et al.| 2008 Jones et al. 


2015). However, in Figure 10 it can be seen that harder 


ionizing spectra give higher values of O 32 at fixed ioniza¬ 
tion parameter. We must then consider the possibility 
that high- and low-metallicity star-forming regions have 
similar ionization parameters while having higher O 32 
values because of the change in the hardness of the ion¬ 


izing spectrum with metallicity. At fixed the range 
of ionizing spectra considered in the models spans ~ 0.9 
dex in O 32 . Local star-forming galaxies from SDSS span 
a range of ~ 1.9 dex in O 32 from log( 032 ) ^ — 1.0 to 0.9 
(see Fig. [^. Based on physically-motivated input spec¬ 
tra, it is not possible for the models to span a range that 
large in O 32 at fixed U. It would require spectra that are 
harder or softer than what is reasonably expected from 
models of the young stellar populations responsible for 
ionizing HII regions. Because the hardness of the ion¬ 
izing spectrum increases with decreasing metallicity and 
low-metallicity objects are observed to have higher O 32 
values on average, the only way to produce the dynamic 
range in O 32 observed among the SDSS sample is for 
high-metallicity objects to have low ionization parame¬ 
ters, and vice versa. The relation between 0/H and ion¬ 
ization parameter must exist in the local universe, along 
with a changing ionizing spectrum with metallicity, to 
reproduce the observed range in O 32 . 

Once established, the existence of an anti-correlation 
between O/H and ionization parameter lends insight into 
the inte rpretation of the simple models presented in Fig¬ 
ure 0 Using si milar photion i zation models to those 
presented above, Steidel et al. (2014) argued that the 
sequence of galaxies in the [IN iij BF i’ diagram at low 
or high redshift is primarily a sequence in ionization 
parameter because increasing ionization parameter with 
all other inputs fixed moves grid points along the star¬ 
forming sequence. Given the dependence of the ioniza¬ 
tion parameter and hardness of the ionizing spectrum 
on 0/H, many of the grid points shown in Figure 11 
are not descriptions of real objects. We do not expect 
to see galaxies with high metallicities, high ionization 
parameters, and hard ionizing spectra. Likewise, galax¬ 
ies with low metallicities, low ionization parameters, and 
soft ionizing spectra are not likely to be observed. In 
fact, it has been previously observed that local HII re¬ 
gions and star-forming galaxies only occupy a narrow 
subset of t he parameter space in photoionization mode l 
grids (e.g., Dopita fc Evaiis|1986 Dopita et al.|2006a|b ). 
While individual HII regions can simultaneously demon- 
strate low metallicity and l ow ionization param e ter, a s 
in the sample discussed in van Zee & Haynes (2006), 
such objects would contribute negligibly to a luniinosity- 
weighted galaxy-averaged spectrum for galaxies that are 
actively star-forming. 

The star-forming sequence in the [N ii] BPT diagram 
can be understood as a sequence in both ionization pa¬ 
rameter and metallicity because the two are fundamen - 


Steidel et al. 


tally linked (Bresoli n et al. 2012 Sanchez et al. 2015). 
^ (2014[) did mention the relations between 


ionization parameter, metallicity, and ionizing spectrum 
as a way to reconcile the utility of strong-line metallicity 
indicators in the local universe with the apparent lack 
of dependence of the star-forming sequence on metallic¬ 
ity in photoionization models. Our results suggest that 
these relations hold in similar form at z ^ 2.3 as well. 
We note that these relations are traced by the typical 
properties of the galaxy population, and do not preclude 
the existence of a small number of galaxies in unexpected 
regions of the parameter space due to scatter in the rela¬ 
tions or significantly different conditions caused by some 
rare process or event, such as a major merger. 
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Figure 12. O32 vs. R23 (left), 03N2 (middle) and N2 (right) for local star-forming galaxies (gray histogram) and the z ~ 2.3 [N ll] 

subsample (red points). In all panels, the white line shows the running median of the local sequence of star-forming galaxies. The dark red 
line shows the running median of the z ~ 2.3 sample in the left panel, and the best-fit linear relation in the middle and right panels. 


5.4. Nitrogen abundance and the cause of the [N n] 
BPT diagram offset 

A significant amount of effort has been put forth to find 
the cause of the offset high-redshift galaxies display in the 


IN III BPT diagram (IKewlev et al.|2013a|b 

Masters et al. 

2014| 

bteidel et ai.||2014 

bhapley et al.|| 2 U] 

lo). Proposed 


eters (Brinchmann et al. 2008 

), systematically harder 

ionizing spectra (|Steidel et al.J 

2014), eleyated nitrogen 

abundance at hxed metallicity (p 

blasters et al. 2014 

Shap- 


(|Kewley et al.|2(J13a). It has additionally been proposed 
that the oftset in the BPT diagram is an artifact arising 
from widesp read, weak AGN in t he high-redshift g alaxy 


population (Wright et al. 2010). Coil et al. (2015) have 
shown that such global ACIN contamination does not ap¬ 
pear to be present among z ~ 2.3 galaxies from the MOS¬ 
DEF survey. 

In this paper, we have sought to characterize the phys¬ 
ical properties influencing the ionization state of high- 
redshift galaxies, and can now use the results presented 
herein to inve stig ate the cause of the [N ii] BPT diagram 


offset. Figure 13 shows the [N il] and [S ii] BPT diagrams 
(Fig.^j) and We O 32 vs. R 23 , 03N2, and N2 diagrams 
(Fig.^^ with the z ^ 2.3 galaxies color-coded according 
to the magnitude of the offset in the [N ii] BPT diagram. 
We divided the z ~ 2.3 sample in the [N ii] BPT diagram 
at the running median [O iil]/H/3 in bins of [N il]/Ha, 
shown by the solid purple line in the top right panel. 
Galaxies above and to the right of this line in the [N ii] 
BPT diagram are shown in blue, while those galaxies 
that show a smaller offset and overlap the local sequence 
are shown in red. We have verified that our results do 
not change if other dividing lines in the [N ii] BPT dia¬ 
gram are used, including a linear fit and t he best-fit line 
to the z ^ 2.3 star-forming sequence from|Shapley et al.| 


(2015) (red dashed line). 

We find that the blue and red data points are well 
mixed and are not systematically offset from the local 
distribution in the O 32 vs. R 23 diagram. The running 
median of the three samples follow a very similar se¬ 
quence where there is overlap at log( 032 ) ~ —0.1 to 0.4. 
Similarly, the blue data points do not appear to follow a 


different distribution from that of the red data points or 
the low-metallicity local sequence in the [S ii] BPT dia¬ 
gram. There is some separation in the two subsamples 
in the [S ii] BPT diagram because the blue data points 
have higher [O iii]/H/? a priori due to the selection. In 
the 03N2 diagram, the red points fall close to the local 
distribution, while the blue points appear to be slightly 
more offset towards lower 03N2 at hxed O 32 . The blue 
data points display less overlap with the red data points 
in the N2 diagram and a larger offset compared to local 
galaxies. We use the best-fit linear relations to the blue 
and red datapoints, respectively, to quantify the mean 
offsets and make the offsets clearer. The red line has a 
mean offset of -0.11 dex in 03N2 at fixed O 32 compared 
to the local median, while the blue line is offset -0.30 dex 
on average. In the N2 diagram, the red and blue lines 
are offset 0.08 and 0.40 dex, respectively, towards higher 
N 2 at fixed O 32 . 

To summarize, we find that, as a whole, the z ^ 2.3 
sample shows no significant offset in line ratio diagrams 
that do not include nitrogen ([S ii] BPT and O 32 vs. 
R 23 diagrams). We find a slight offset in a diagram that 
utilizes a line ratio including nitrogen as well as oxygen 
(O 32 vs. 03N2 diagram). This offset increases in the N2 
diagram, in which nitrogen is the only metal species in 
the line ratio. Additionally, we hnd that galaxies showing 
the largest offset in the [N il] BPT diagram do not distin¬ 
guish themselves from galaxies showing a smaller offset 
when plotted in diagrams that do not include nitrogen. 
The most-offset galaxies in the [N il] BPT diagram dis¬ 
play systematically lower 03N2 values at fixed O 32 than 
galaxies showing a small offset in the [N il] BPT diagram. 
This effect increases in magnitude in the O 32 vs. N2 di¬ 
agram, in which those galaxies showing a large offset in 
the [N il] BPT diagram display signihcantly larger N2 
values at fixed O 32 than the less-offset subsample of the 
z ^ 2.3 galaxies. A simple way to change the [N ii]A6584 
flux without affecting the flux of the other strong optical 
emission lines is to change the nitrogen abundance. Our 
observations suggest that the N/0 ratio at fixed 0/H is 
higher on average in z ~ 2.3 star-forming galaxies com¬ 
pared to local star-forming galaxies. The high-redshift 
galaxies showing the largest offset in the [N il] BPT di- 
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Figure 13. Emission line ratio diagrams with z ~ 2.3 galaxies color-coded by the offset in the [N ll] BPT diagram. In each panel, the 
gray histogram shows the distribution of local star-forming galaxies from SDSS. The top row shows the [N ll] (top left) and [S ll] (top right) 
BPT diagrams. The O 32 vs. R 23 , 03N2, and N2 diagrams are shown in the left, middle, and right panels, respectively, of the bottom row. 
In each panel, the white line shows the running median of the 2 : ~ 0 sample. In the top left panel, the solid purple line shows the running 
median of the z ~ 2 .3 star-forming galaxies, while the dashed red line shows the fit to the MOSDEF z ~ 2.3 star-forming sequence from 
IShapley et al.| ||2015|l. In all panels, z ~ 2.3 galaxies falling above and to the right of the dashed red line are color-coded blue, while z ~ 2.3 
galaxies tailing below and to the left of this line are color-coded red. In the lower left panel, the red and blue curves show the median R 23 
in bins of O 32 for the red and blue data points, respectively. The solid red and blue lines in the lower middle and right panels show the 
best linear fit to the red and blue data points, respectively. See Footnote 9 regarding the use of running medians and best-fit lines. 


agram appear to have higher N/0 ratios than less offset 
z ~ 2.3 galaxi es of the s ame metallicity. 

In Sections |5.1| and |5.2[ we showed that our data 
are inconsistent with a systematically harder ionizing 
spectrum or higher ionization parameter in high-redshift 
galaxies at fixed metallicity. An increase in the gas den¬ 
sity or ISM pressure can also move galaxies in the di¬ 
rection of the observed offset. We utilize Cloudy pho¬ 
toionization models to quantify the magnitude of this 
effect using the characteristic densities obtained in Sec¬ 
tion Figure [T4| shows the local star-forming sequence 
for rraerence, along with two photoionization model grids 
produced with the same ionizing spectrum (blackbody 
with Teff = 50, 000 K) while varying the electron density 
between 25 cm“^ and 250 cm“^. The grids are shown 
over the metallicity range 0.2-0.6 Zq in which the typi¬ 
cal metallicity at z ~ 2.3 is expected to fall. The low¬ 
est metallicity grid points show negligible change in line 
ratios with density, while the dependence on density in¬ 
creases with metallicity. With an increase of a factor of 
10 in density, [O iii]/H/3 and [N Ii]/Ha are increased by 


<0.1 dex at fixed ionization parameter and metallicity. 
This is not a large enough shift to account for the ob¬ 
served offset in the [N il] BPT diagram. If our density 
estimate for local star-forming regions is underestimated, 
then the magnitude of the line ratio shift caused by den¬ 
sity will be even smaller. Thus, the increase in density 
from z^0toz^2.3 plays only a minor role in the [N ii] 
BPT diagram offset. 

We conclude that the observed offset of z ^ 2.3 galax¬ 
ies from the local star-forming sequence in the [N il] BPT 
diagram is mostly caused by elevated N/0 at fixed O/H 
compared to local galaxies, while an increase in den¬ 
sity/pressure of the star-forming regions plays a minor 
role. We do not find evidence that a change in the ioniz¬ 
ing spectrum or ionization parameter at fixed metallicity 
plays a part in the offset. This finding is cons i stent w ith 
earlier MOSDEF results from Shapley et al. (20151 at¬ 
tributing the offset in the BP'i' diagram to higher N/0 
in high-redshift galaxies with M*< 10^° Mq. O ur results 
are also in ag reement with the interpretation of Masters) 


et al. (2014), in which anamolous nitrogen abundance 
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- 1.2 - 1.0 - 0.8 

log([N ii]/Ha) 

Figure 14. Model grids assuming a gas density of 25 cm“^ (light 
green) and 250 cm~^ (cyan) with gas-phase metallicity spanning 
0.2-0 .6 Z 0 and an input ionizing spectr um of a 50,000 K blackbody. 
Grid points are displayed as in Figure [TT] 


was fir st proposed as a cause of the offset. [Steidel et al.| 
(|2014|) also found an elevated N/0 among their z ~ 2.3 
sample, but argued that the primary cause of the off¬ 
set lies in harder ionizing spectra and higher ionization 
parameters in comparison to what is observed locally. 

The cause of the high er N/0 ratios observ ed at high- 
redshift is not yet clear. Masters et al. (2014) speculated 
that a larger-than-normal population of Wolt-Rayet stars 
could produce a nitrogen enhancement, although a mech¬ 
anism to produce such a population exclusively at high 
redshifts was not proposed. If the nitrogen enhancement 
is due to a significant difference in the stellar population 
at a given metallicity, there would likely be a difference in 
the overall ionizing spectrum that population produces, 
for which we do not see evidence. Another possibility 
is that gas flows lead to a larger nitrogen abundance at 
fixed 0/H. If a large amount of unenriched gas is ac¬ 
creted and mixed into the star-forming regions before 
intermediate-mass stars enrich the ISM with nitrogen, 
O/H will decrease. The intermediate-mass stars will then 
release an amount of nitrogen that is larger than what 
is expected at that 0/H because the gas is now less en¬ 
riched than the stars. This scenario, or another mecha¬ 
nism involving gas inflows and outflows, is perhaps more 
likely in high-redshift galaxies which are known to have 
large gas fractions and concentrated star formation that 
could correspond to active accretion of metal-poor gas 
( Tacconi et al.|2(JT3 ). We found evidence consistent with 
an accumulation of unenriched gas in z ^ 2.3 galaxies 


throug h the mass-metallicity relation in Sanders et al. 
(2015). Observations of the rest-frame UV could de- 
termme if exotic stellar populations are present, while 
more observations of the cold gas content of high-redshift 
galaxies with, e.g., ALMA could help uncover the role of 
gas flows. 

The presence of AGN within the z ~ 2.3 sample could 
significantly bias conclusions drawn from excitation di¬ 
agrams. We did not use a locally-determined demarca¬ 
tion to separate star-for ming galaxies from AGN in the 
[N ii] BPT diagram (e.g., Kewley et al.|2001||Kauffmann 
et al. 2003) because such a selection would introduce a 
bias against objects with high N/0 at fixed [O iil]/H/3. 
While some of the galaxies in the z ~ 2.3 sample fall 


in the local AGN region of the [N ii] BPT diagram, we 
do not see any z ~ 2.3 galaxies present where the AGN 
sequence would fall in the other four panels of Figure [T^ 
Therefore, we are confident that our z ~ 2.3 sample is 
free of AGN. 

5.5. Diffuse ionized gas and the interpretation of global 
galaxy speetra 

The observed positions of local and z ~ 2.3 galaxies 
in the [N il] and [S ii] BPT diagrams are key pieces of 
evidence in our arguments regarding the hardness of the 
ionizing spectrum of loca l and high-r eds hift star-forming 
regions (see Section 5.1 and Figure 11). One source of 
uncertainty in the interpretation of these plots is the in¬ 
clusion of light from the diffuse ionized component of the 
ISM in integrated-light galaxy spectra. Emission from 
diffuse ionized gas can be significant, with roughly half 
of the total Ha emission of local spiral galax ies coming 
from a diffuse component ( Zurita et al.|200() ). If the dif¬ 
fuse ionized component has line ratios that are not equiv¬ 
alent to those of HII regions, which is probable because 
diffuse gas and HII regions are characterized by different 
ionizing sources and ionization states, then light from 
the diffuse component can act as a contaminant when 
investigating properties of star-forming regions. This is¬ 
sue is especially concerning for low ionization states such 
as [N ii] and [S ii] which are easily ionized and can have 
elevated flux compared to that o f Ha when shock-excited 
( Martin|[l997[ Hong et al.||2013 ). Gharacterizing the im- 
pact of diffuse emission requires knowledge of both the 
fraction of emission and the emission line ratios coming 
from the diffuse ionized component, at both low and high 
redshift. 

In the local universe, comparison of integrated galaxy 
spectra with those of individual HII regions can give in¬ 
sight into the effect of emission from other components 
of the ISM. Local HII regions have been found to follow 
a similar sequence to that of local star-forming galax¬ 
ies in the [N ii] BPT diagram, but appear to have sys¬ 
tematically lower [S ii]/Ha at fixed [O Iiij/H/I in the 


[S ii] BPT diagram dVeilleux fc Osterbrock||1987 


gin et al. 2012[ Berg et al. ' 21JI5p Groxall et al. 2(JI5p . 


Pilyu- 


Every emiission line in an integrated galaxy spectrum 
contains a luminosity-weighted contribution from each 
HII and diffuse emission region falling in the slit or fiber, 
complicating the comparison of HII region and integrated 
galaxy spectra. Ongoing spatially-resolved sp ectroscopic 
surve ys, such as the SDSS-IV/MaN GA IFU (Law et al. 
20151 and SAMI (Bryant et al.|2015) galaxy surveys, will 
provide a dataset capable of unraveling the relative con¬ 
tributions of HII regions and diffuse ionized gas to global 
galaxy spectra. 

The situation is much more uncertain at high redshifts 
because the structure of the ISM in z ~ 2 galaxies is 
poorly constrained. The ISM is almost certainly differ- 
ent in z ~ 2 galaxies, considering the high gas fractions 
dTacconi et al. 2013) and compact sizes (van der Wei 
et al.||2014|), combi ned with high rates of star formation 
dShivaei et al.|2Ul^, significant feedback driving outflows 
(Shapley et al.||2003| |Steidel et al.||2010P, and hi gh level 
of turbulence observ ed m z ^ 2 disk galaxies (Forsterj 
Schreiber et al.|2009 ). If high-redshift galaxies are domi¬ 
nated by giant kpc-scale HII regions ionized by super-star 
clusters that fill a large fraction of the galaxy volume. 
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the filling factor of the diffuse ionized component may 
be small along with the fraction of line emission origi¬ 
nating there. In this case, it may be more appropriate 
to compare the integrated emission-line spectra of high- 
redshift galaxies to those of luminous local HII regions 
rather than global galaxy spectra. We are only starting 


structure of the ionized ISM at z ~ 2 (Forster Schreiber 

et al .|20111 

Genzel et al.|20111|Jones et al.|2013t JNewman 

et al.| 2014 

). Adaptive optics observations of increased 


the spatially-resolved structure of the ionized component 
of the ISM at high redshift, including the strength of the 
[N il] and [S ii] emission lines. 


6. SUMMARY 

We have investigated the physical conditions of star¬ 
forming regions at z ~ 2.3, specifically the electron den¬ 
sity and ionization parameter, and made comparisons to 
local star-forming galaxies in order to understand how 
these properties evolve with redshift. We performed this 
investigation using rest-frame optical spectra of z ~ 2.3 
galaxies from the ongoing MOSDEF survey. We sum¬ 
marize our main conclusions below and discuss future 
observations that could shed additional light on the ion¬ 
ization state of high-redshift galaxies. 


1. We explored the evolution of the electron density 
in star-forming regions using the [O ii]AA3726,3729 
and [S ii]AA6716,6731 doublets. We found that 
z ^ 2.3 galaxies have median [O ii]A3729/A3726 = 
1.18 and median [S ii]A6716/A6718 = 1.13, cor¬ 
responding to electron densities of 225 cm“^ and 
290 cm“^, respectively. Local star-forming galaxies 
from SDSS have median [S ii]A6716/A6718 = 1.41 
which yields a density of 26 cm“^. We found an 
evolution in electron density of an order of magni¬ 
tude between z ~ 0 and z ~ 2.3. 


2. We investigated the ionization state of z ^ 2.3 and 
local star-forming galaxies, using O 32 as a proxy for 
the ionization parameter. We found that O 32 de¬ 
creases with increasing stellar mass in both the lo¬ 
cal and z ^ 2.3 samples. The slope of the O 32 —M* 
anti-correlation is very similar for both samples, 
but the z ~ 2.3 sample is offset ^0.6 dex higher 
in O 32 relative to the local sample at fixed M*. 
This offset can be explained by the evolution of 
the mass-metallicity relation with redshift, such 
that high-redshift galaxies have lower metallicities 
at fixed M*, and the existence of an anti-correlation 
between ionization parameter and 0/H. 

3. We found that z ^ 2.3 galaxies show no systematic 
offset from local galaxies and follow the distribu¬ 
tion of the low-metallicity tail of local galaxies in 
the O 32 vs. R 23 diagram. The high-redshift sam¬ 
ple behaves similarly to the local sample in the O 32 
vs. 03N2 diagram, displaying a slight systematic 
offset from the local distribution. We propose that 
z ~ 2.3 galaxies follow the same anti-correlation 
between ionization parameter and 0/H that is ob¬ 
served in the local universe. 


4. Using simple photoionization models, we demon¬ 
strated that the translation between O 32 and ion¬ 
ization parameter is only strongly dependent on the 
shape of the ionizing spectrum, and has little de¬ 
pendence on the assumed gas density and gas-phase 
metallicity. This translation will only evolve with 
redshift if the hardness of the ionizing spectrum at 
fixed metallicity evolves with redshift. 

5. We utilized the position of z ^ 2.3 and local galax¬ 
ies in the [N ii] and [S ii] BPT diagrams combined 
with simple photoionization models to show that 
the hardness of the ionizing spectrum does not sig¬ 
nificantly increase or decrease with redshift. Pho¬ 
toionization models predict that a hardening of the 
ionizing spectrum will increase [O iii]/H/3 at fixed 
[N ii]/Hq! and [S il]/Ha. The z ^ 2.3 sample dis¬ 
plays an offset from the local sequence in the [N il] 
BPT diagram but is not significantly offset in the 
[S il] BPT diagram. We conclude that there is not 
a significant increase in the hardness of the ioniz¬ 
ing spectrum at fixed metallicity between z ^ 0 
and z ~ 2.3. 

6 . Galaxies at z ^ 2.3 show no significant systematic 
offset from local galaxies in line ratio diagrams in¬ 
volving only lines of oxygen, sulfur, and hydrogen, 
while they show a systematic offset in line ratio di¬ 
agrams involving nitrogen. These results suggest 
that metallicity indicators using line ratios exclud¬ 
ing nitrogen (e.g. the combination of O 32 and R 23 ) 
do not evolve up to z ^ 2.3, while indicators using 
nitrogen are biased due to an evolution in N/0 at 
fixed 0/H. 

7. A consequence of conclusions 6 and 7 is that z ~ 2.3 
have similar ionization parameters to z ~ 0 galax¬ 
ies at fixed metallicity. Higher typical ionization 
parameters are inferred for z ~ 2.3 galaxies com¬ 
pared to those of local galaxies because z ^ 2.3 
galaxies have lower typical metallicities. The ion¬ 
ization state appears to be set by the metallicity 
both locally and at z ~ 2.3. 

8 . We investigated the offset between z ^ 2.3 and lo¬ 
cal galaxies in the [N il] BPT diagram. We found 
that the z ~ 2.3 galaxies that display the largest 
offsets in the [N il] BPT diagram are not signifi¬ 
cantly offset from the local distribution or the re¬ 
mainder of the z ^ 2.3 sample in the O 32 vs. R 23 
and [S ii] BPT diagrams, but are systematically 
offset in the O 32 vs. 03N2 and N2 diagrams. We 
conclude that higher N/0 at fixed 0/H drives the 
z ^ 2.3 offset in the [N ii] BPT diagram. We pre¬ 
viously provided evidence against significant evo¬ 
lution of the hardness of the ionizing spectrum or 
ionization parameter at fixed 0/H. We used simple 
photoionization models to show that an evolution 
of a factor of 10 in the gas density cannot account 
for the full offset in the [N ii] BPT diagram. We 
further conclude that an increase in the gas density 
plays a minor secondary role in driving the z ~ 2.3 
offset in the [N ii] BPT diagram. 

There still remain many questions to be answered re¬ 
garding the ionization state of local and z ^ 2.3 galax- 
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ies. Emission line contribution from the diffuse ion¬ 
ized component of the ISM is uncertain both locally 
and at high redshifts. Variety in photoionization mod¬ 
eling leads to variation in estimated ionization param¬ 
eters. There are some observations that would allow 
us to test our proposed scenario that z ^ 2.3 galaxies 
have the same ionization parameter as z ^ 0 galaxies 
with the same metallicity. The most obvious of these 
is direct-method oxygen abundance measurements from 
auroral lines at z ^ 2. Direct-method abundances are 
insensitive to changes in the electron density, showing 
variations of < O.OI dex from 25 cm“^ to 250 cm“^ 
(i.e., the observed evolution in density from z ^ 0 to 
z ~ 2.3), an advantage over some strong-line methods 


(see Jones et al. 20151. Currently, only a handful of 
direct-method abundance measurements have been at- 


(Villar-Martm et al. 2004[ 

Yuan & Kewley 2009 [Erb 

et al.||2010l IRigby et al.| 

120111 IHrammer et ai. 2012bt 

Christensen et ai. 2012 |. 

Vlaseda et al. 2014). Current 


scopes and upcoming instruments on the Thirty Meter 
Telescope will allow for observations of auroral lines for 
more typical galaxies at z > 1. Additionally, obser¬ 
vations of other ionization-parameter-sensitive emission 
line ratios could provide a test o f this scenario, including 
[Ne ii i]A3869/[0 il]AA3726,3729 ([Levesque fc Richardson 


2014| and [S iii ]AA9069,9532/[S iiJAA6716,6731 (Kewley 


&: Dopita]|2002 ). Finally, improved stellar and photoion- 


ization models that can make self-consistent predictions 
of the physical properties of star-forming galaxies across 
the entire local star-forming sequence are needed to cre¬ 
ate more reliable ionization parameter calibrations. Such 
models would provide a local foundation for reliable pho¬ 
toionization models of the full suite of strong rest-frame 
optical emission lines out to z ~ 2. 
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